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ABSTRACT 


This  thesis  was  motivated  by  the  possibility  of  using  high  energy  laser  pulses  as 
an  acoustic  source  for  naval  applications.  Research  conducted  in  the  1970’s  and  80’s 
shows  that  sound  production  from  laser  pulses  is  most  efficient  when  the  energy  density 
of  the  pulse  exceeds  the  threshold  required  for  plasma  formation.  The  resulting  acoustic 
wave  falls  into  the  highly  non-linear  shock  regime.  Later  work  by  Vogel  et  al.  sought  a 
more  complete  understanding  of  the  non-linear  dynamics  and  energy  distribution  of  this 
process  in  an  attempt  to  limit  collateral  tissue  damage  during  laser  surgery.  This  work 
includes  detailed  experimental  data  including  photographs  and  hydrophone 
measurements  as  well  as  numerical  calculations  of  expected  pressures,  bubble  dynamics, 
and  pulse  shapes. 

The  goal  of  this  thesis  was  to  investigate  the  characteristics  of  the  laser  generated 
acoustic  pulse  further  through  experimentation  and  modeling.  Experiments  were  carried 
out  with  Ted  Jones  at  the  Naval  Research  Laboratory  to  investigate  the  directionality  of 
the  acoustic  pulse  produced  by  a  lOOfs  2mJ  laser  pulse  focused  just  under  the  surface  of 
water.  The  range  dependence  of  the  pressure  amplitude  was  also  examined.  The 
amplitude  of  the  pulse  was  found  to  vary  with  direction;  however,  this  effect  is 
considered  likely  to  be  a  result  of  interference  between  the  direct  path  and  the  surface 
reflection.  A  linear  least-squares  fit  of  the  peak  pressure  amplitude  with  range  revealed  a 
1/r1'2  relationship.  This  is  consistent  with  the  expected  approximately  1/r  relationship  for 
pressure  amplitudes  under  lOOMPa.  The  modeling  effort  employed  AUTODYN,  a  finite 
element  program  designed  to  handle  the  non-linear  processes  in  explosions.  The  laser 
generated  acoustic  source  was  modeled  using  an  explosive  of  the  same  volume  as  the 
laser  spot  reported  by  Vogel  for  his  lOmJ  6ns  pulse.  The  internal  energy  of  the  explosive 
was  adjusted  until  the  pressure  amplitudes  agreed  with  Vogel’s  measured  values.  The 
efficiency,  pulse  length,  pulse  shape,  and  variation  of  pressure  amplitude  with  range 
achieved  with  AUTODYN  are  comparable  to  those  reported  by  Vogel. 
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I.  INTRODUCTION 


A.  MOTIVATION 

One  of  the  biggest  threats  to  a  naval  ship  is  a  torpedo  attack  from  an  enemy 
submarine.  Although  an  anti-torpedo  torpedo  (ATT)  is  in  development,  the  only 
protection  ships  currently  have  against  torpedoes  is  a  combination  of  maneuvering  the 
ship  and  the  use  of  acoustic  countermeasures.  Even  after  the  ATT  is  fielded,  acoustic 
countermeasures  will  continue  to  be  an  important  tool  for  surviving  enemy  torpedo 
attacks.  These  countenneasures  can  be  deployed  as  expendable  sources  thrown  over  the 
side  or  as  a  towed  source.  The  AN/SLQ-25  system,  also  known  as  Nixie,  is  an 
underwater  acoustic  source  towed  behind  a  ship  and  connected  to  an  onboard  signal 
generator  by  a  coaxial  cable.  Nixie  is  capable  of  operating  in  a  broadband  mode  to 
simulate  the  passive  signature  of  a  ship.  It  can  also  detect  the  torpedo’s  ping  frequency, 
amplify  it  2  to  3  times,  and  then  retransmit  it  in  an  attempt  to  simulate  the  ship’s  echo  and 
lure  the  torpedo  away  from  the  ship.  The  problem  with  Nixie  is  that  it  takes  at  least  15 
minutes  to  deploy  fully.  If  the  ship  is  in  known  enemy  submarine  territory  this  is  not  a 
problem  because  Nixie  is  most  likely  already  deployed,  but  Nixie  provides  no  protection 
from  a  surprise  attack.  This  suggests  the  advantage  of  having  an  acoustic  countermeasure 
that  is  quick  to  deploy. 

Much  research  has  been  done  over  the  years  on  the  use  of  high  energy  lasers 

(HELs)  as  a  tool  for  anti-missile  defense.  Such  lasers  are  being  investigated  as  a  possible 

replacement  for  the  MK15  Phalanx  Close-In  Weapons  System  (CIWS),  a  shipboard 

missile  defense  system.  In  addition  to  causing  ablation  on  the  surface  of  an  incoming 

missile,  the  energy  in  a  laser  pulse  is  also  capable  of  producing  sound  through  its 

interaction  with  water.  Assuming  that  the  laser  was  ready  to  fire,  the  speed  with  which 

sound  could  be  produced  by  the  laser  far  from  the  ship  would  be  practically 

instantaneous.  Therefore,  it  offers  the  possibility  of  being  a  valuable  tool  for  anti-torpedo 

defense.  Not  only  would  it  be  faster  than  a  towed  system,  but  it  might  also  eliminate  the 

requirement  for  expendable  countermeasures.  To  be  a  feasible  countermeasure  the  laser- 
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generated  source  would  need  to  be  able  to  produce  multiple  frequencies  in  the  kilohertz 
range  and  Sound  Pressure  Levels  in  the  hundreds  of  decibels.  The  purpose  of  this  thesis 
is  to  examine  the  sound  produced  by  a  high  energy  laser  pulse  to  gain  a  better 
understanding  of  whether  it  might  serve  as  a  feasible  acoustic  source  for  torpedo 
countermeasures  and/or  other  shipboard  applications  of  underwater  sound. 

B.  HISTORY 

The  idea  that  a  light  source  could  be  used  to  produce  sound  was  first  proposed  by 
Alexander  Graham  Bell.  In  his  experiments,  he  showed  that  it  was  possible  to  generate 
sound  waves  in  a  medium  that  could  absorb  the  light  energy.  This  process  became  more 
interesting  with  the  invention  of  the  laser.  In  the  1980’s,  research  into  the  use  of  air- 
mounted  high  energy  pulsed  lasers  for  missile  defense  led  naturally  to  the  investigation  of 
this  technology  for  a  dual  use  in  underwater  acoustics.  Some  work  also  appears  to  have 
been  motivated  by  Soviet  research  into  laser-generated  acoustics.  In  any  case,  several 
U.S.  researchers  including  Bruce  Maccabee  from  the  Naval  Surface  Warfare  Center 
Dahlgren,  Ilene  Busch-Visniac  and  Yves  Berthelot  from  the  Applied  Research 
Laboratory  in  Austin,  and  Allan  Pierce  and  H.  Hsieh  from  the  Georgia  Institute  of 
Technology  all  looked  at  the  feasibility  of  using  high  energy  pulsed  lasers  as  acoustic 
sources.  Interest  in  laser-generated  sound  was  also  motivated  in  the  mid-1990s  as  lasers 
became  a  popular  surgery  technique.  This  led  to  experiments  by  A.  Vogel,  S.  Busch,  and 
U.  Parlitz  in  which  they  investigated  the  shock  waves  produced  by  fairly  low  energy  laser 
pulses  in  water.  In  the  late  1990’s,  there  is  also  evidence  of  interest  in  laser-generated 
acoustics  for  naval  purposes  among  the  Chinese.  He  and  Feng  from  the  Chinese 
Academy  of  Sciences  in  Shanghai  and  Duo-Min  from  Nanjing  University  described 
experimental  work  in  which  they  produced  sound  in  water  using  laser  pulses.  Duo-Min 
specifically  mentions  underwater  acoustic  communications  as  a  possible  application. 
Most  recently,  Ted  Jones  from  the  Naval  Research  Laboratory  and  his  collaborators  from 
a  wide  variety  of  government  labs  as  well  as  private  industry  have  developed  some 
interesting  new  techniques  to  increase  the  achievable  sound  levels. 
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c. 


CONTRIBUTIONS  OF  THIS  THESIS 


Experiments  were  conducted  for  this  thesis  at  the  Naval  Research  Laboratory 
along  with  Ted  Jones  and  his  group  at  Code  6795.  These  experiments  involved  the  use  of 
a  titanium-sapphire  laser  producing  femtosecond  pulses  with  an  energy  of  about  2mJ  and 
a  wavelength  of  about  400nm.  An  analysis  of  these  experiments  includes  an  analysis  of 
the  apparent  directionality  of  the  laser-generated  acoustic  pulse  and  its  range  dependence. 

The  next  part  of  this  thesis  involves  an  attempt  to  model  the  laser-generated 
acoustic  pulse.  The  finite  element  program  COMSOL  Multiphysics  was  used  at  first; 
however,  it  is  not  designed  to  handle  non-linear  processes.  For  this  reason,  the  decision 
was  made  to  switch  to  AUTODYN.  AUTODYN  is  also  a  finite  element  program,  but  it 
is  specifically  designed  to  handle  the  non-linear  effects  involved  with  high  energy 
explosions.  When  the  energy  of  a  laser  pulse  is  deposited  quickly  enough  into  a  small 
volume  of  water,  it  evaporates  and  ionizes  the  water.  This  causes  a  rapid  expansion  in 
the  volume  of  the  bubble,  which  produces  a  shock  wave.  AUTODYN  was  used  in  an 
attempt  to  model  this  shock  wave. 

D.  ORGANIZATION 

Chapter  II  of  this  thesis  includes  a  more  detailed  review  of  the  existing 
background  work  on  laser-generated  acoustic  sources.  The  focus  of  this  chapter  is  on  the 
physics  of  this  process  and  the  acoustic  pressures,  which  have  been  achieved.  Chapter  III 
on  Theoretical  Acoustic  Pulse  Estimates  provides  a  qualitative  summary  of  the 
characteristics  of  both  the  laser  pulse  and  the  acoustic  medium,  which  affect  the 
amplitude,  shape,  length,  and  directionality  of  the  acoustic  pulse.  It  also  presents  some 
issues  with  previous  attempts  to  quantify  where  the  initial  electro-optic  energy  ends  up. 
Chapter  IV  follows  with  a  description  of  the  experimental  set-up  for  the  work  done  at  the 
Naval  Research  Laboratory  and  an  analysis  of  the  results.  Chapter  V  explains  the 
modeling  effort.  It  goes  through  the  various  approaches  tried  in  AUTODYN  to  model 
the  shock  wave  generated  by  the  laser  pulse.  This  discussion  also  includes  some  insights 
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into  the  limitations  of  AUTODYN  for  this  particular  application.  Chapter  VI  discusses 
the  conclusions  to  be  drawn  from  this  work,  and  recommendations  for  future  work  are 
provided  in  Chapter  VII. 
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II.  BACKGROUND 


In  this  chapter,  the  background  work  mentioned  in  the  introduction  will  be 
explored  in  greater  detail.  An  overview  of  the  relevant  work  will  be  provided  as  well  as 
details  concerning  the  theory  developed  and  the  experimental  values  obtained. 

A.  OVERVIEW  OF  PREVIOUS  WORK 

During  the  1980’s,  there  was  considerable  interest  in  high  energy  lasers  as  a 
possible  technology  for  missile  defense.  In  particular,  the  feasibility  of  mounting  high 
energy  lasers  on  board  aircraft  was  examined.  This  interest  raised  the  question  as  to 
whether  the  lasers  could  be  used  for  other  applications  as  well.  Building  on  the 
theoretical  and  experimental  work  of  both  US  and  Soviet  research  in  the  late  1970’s,  at 
least  three  U.S.  groups  looked  at  this  issue.  Yves  Berthelot  and  Ilene  Busch-Visniac 
from  the  Applied  Research  Laboratories,  The  University  of  Texas  at  Austin  [1985] 
investigated  the  acoustic  wave  generated  by  a  laser  source  moving  over  the  water  at 
speeds  ranging  from  subsonic  and  supersonic.  Allan  Pierce  and  H.  A.  Hsieh  of  the 
Georgia  Institute  of  Technology  [1986]  proposed  an  approach  designed  to  tailor  the 
directionality  of  the  acoustic  signal  by  the  careful  placement  and  timing  of  laser  beams  on 
the  ocean  surface.  Bruce  Maccabee  [1987]  from  the  Naval  Surface  Warfare  Center 
Carderock  Division  provided  the  theory  for  linear  pressure  waves  resulting  from  simple, 
single  laser  pulses  as  well  as  experimental  data  on  the  pressures  obtained.  The  interest  in 
airborne  laser  systems  diminished  in  the  1990’s  most  likely  due  to  a  drop  in  the  funding 
for  this  form  of  strategic  missile  defense;  however,  interest  in  laser  generated  sound 
resurfaced  in  the  mid  to  late  1990’s  with  the  increasing  popularity  in  laser  surgery. 
German  scientists  A.  Vogel,  S.  Busch,  and  U.  Parlitz  [1996]  published  an  article  with  a 
theoretical  treatment  of  the  shock  waves  produced  by  the  laser  pulses  in  water  as  well  as 
their  experimental  observations.  This  work  was  motivated  by  a  desire  to  understand  the 
potential  tissue  damage  attendant  to  laser  surgery  -  particularly  in  the  eye.  Chinese 
researcher,  He  Duo-Min  of  the  Nanjing  University  of  Science  and  Technology  showed 
renewed  interest  in  the  use  of  laser  generated  sound  for  submarine  communications  in  a 
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1998  paper  delivered  at  a  Strategic  Research  Seminar  in  Singapore.  Two  other  Chinese 
researchers,  H.  He  and  S.  Feng  at  the  Chinese  Academy  of  Sciences  in  Shanghai  [1999], 
also  showed  an  interest  in  laser  generated  acoustics;  however,  they  were  using  the 
phenomenon  as  a  means  of  investigating  the  optical  response  time  of  the  water  through 
its  ability  to  focus  an  intense  laser  beam  by  changes  in  its  refractive  index.  Also  in  the 
late  1990’s,  Ted  Jones,  Jacob  Grun,  Ray  Burris,  and  Charles  Manka  [1999]  published  a 
technical  report  at  the  Naval  Research  Laboratory  (NRL)  which  was  to  be  the  starting 
point  for  a  program  designed  to  explore  and  improve  the  acoustic  signal  which  could  be 
generated  using  high  energy  laser  pulses.  By  2006,  this  effort  led  to  a  novel  idea  for 
increasing  the  intensity  of  the  laser  pulse  in  water,  i.e.,  the  use  of  a  reverse  chirp  in 
conjunction  with  dispersion  to  shorten  the  pulse  in  the  direction  of  propagation. 

Of  these  early  works,  several  will  be  explored  in  more  detail.  The  basic  theory  of 
the  laser- water  interaction  will  be  presented  for  both  linear  and  shock  waves.  In  addition, 
the  basic  ideas  for  improving  the  performance  of  a  laser  acoustic  source  will  be  explained 
and  a  comparison  of  the  experimental  data  obtained  by  various  groups  will  be  provided. 

B.  MECHANISMS  FOR  LASER-GENERATED  SOUND 

Maccabee  [1987]  proposed  that  the  mechanisms  for  laser  generated  sound  be 
divided  into  four  types.  It  is  important  to  distinguish  these  different  types,  because  the 
theory  as  well  as  the  efficiency  of  the  process  varies  substantially  depending  on  the 
nature  of  the  laser-water  interaction.  The  first  type  occurs  when  the  laser  energy  heats 
the  water  but  is  not  high  enough  to  cause  evaporation.  The  electro-optic  to  acoustic 
efficiency  for  this  type  of  interaction  is  approximately  0.01%.  The  second  type  occurs 
when  the  absorbed  energy  is  high  enough  to  cause  evaporation.  This  rapid  evaporation 
causes  the  volume  to  expand  very  fast  and  produces  a  larger  pressure  than  in  the  first 
type.  The  efficiency  for  Type  2  is  higher  than  the  Type  1  at  0.1-1%.  The  third  type  is 
where  the  laser  energy  is  so  high  that  it  ionizes  the  water.  The  formation  of  a  plasma 
allows  the  water  to  absorb  more  of  the  laser  energy,  thus  causing  the  volume  to  expand 
faster  and  to  a  larger  final  radius.  This  type  has  the  highest  efficiency  at  10-30%.  The 
final  type  of  laser  generated  acoustics  is  produced  by  focusing  the  laser  onto  an 
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underwater  object.  For  this  case,  the  frequency  of  the  laser  should  be  one  that  has  the 
lowest  absorption  coefficient  in  water  so  that  most  of  the  laser  energy  will  be  absorbed  by 
the  underwater  object  instead  of  the  water.  From  here,  the  process  is  similar  to  type  one. 
The  absorption  of  laser  energy  will  cause  the  object’s  temperature  to  increase,  which  will 
cause  it  to  expand  rapidly  producing  a  pressure  wave.  Since  the  efficiency  of  the  Type  4 
process  depends  on  the  particulars  of  the  object,  Maccabee  does  not  speculate  on  the 
range  of  efficiencies  which  might  be  expected  from  it. 

Maccabee  provides  some  interesting  experimental  data  of  the  pressure  waves 
produced  by  laser  pulses.  This  data  was  taken  4.9  cm  below  the  laser  pulse,  and  the 
results  are  shown  in  the  graph  below.  The  cross-sectional  area  of  the  spot  remained 
constant  at  0.2cnr.  The  pulse  length  was  changed  between  the  third  and  fourth  graphs, 
and  measurements  were  not  made  of  the  spot  volume.  Thus,  the  exact  energy  densities 
are  unknown.  However,  it  is  clear  from  the  plots  that  at  low  laser  pulse  energies,  the 
acoustic  pulse  assumes  a  bipolar  shape  while  at  higher  energies  it  starts  to  show  a  shock 
front.  The  acoustic  pressure  goes  below  an  absolute  pressure  of  zero.  This  is  possible 
only  when  the  timeframe  of  the  process  is  too  quick  for  evaporation  to  occur,  i.e.,  the 
water  is  in  tension.  At  yet  higher  pulse  energies,  the  negative  pressure  swing  disappears 
and  a  shock  wave  with  exponential  tail  appears. 
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Figure  1  Experimental  results  of  the  pressure  wave  received  4.9  cm  below  laser  spot 
for  varying  amounts  of  laser  pulse  energies.  (From:  Maccabee,  1987) 


Of  the  types  of  laser  sound  production  mentioned  by  Maccabee,  Type  1  (heating) 
and  Type  3  (ionizing)  are  of  the  most  interest.  Type  1  is  of  interest  because  it  can  be 
described  by  linear  theory.  Type  3  is  of  interest,  because  the  efficiency  of  the  electro¬ 
optic  to  acoustic  conversion  is  substantially  higher  than  Type  1  or  2.  Type  4  (focusing 
energy  onto  an  object)  is  potentially  interesting,  but  the  details  would  depend  on  the 
nature  of  the  object  used.  Therefore,  only  Type  1  and  3  interactions  will  be  examined  in 
more  detail. 

C.  TYPE  1  SOUND  PRODUCTION  -  LINEAR  THEORY 

When  a  laser  pulse  deposits  energy  in  water  such  that  the  energy  density  is  below 
the  threshold  required  for  evaporation,  the  water  expands  proportionally  to  the 
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temperature  increase  and  the  volume  coefficient  of  expansion.  The  theory  for 
determining  the  pressure  and  directionality  of  the  resulting  pressure  wave  is  identical  for 
any  mechanism  which  heats  water.  Sulak  et  al.  [1979],  for  example,  derive  an  expression 
for  the  acoustic  wave  resulting  from  the  interaction  of  a  proton  beam  with  water. 
Interestingly,  they  point  out  that  at  temperatures  below  6°C,  the  coefficient  of  expansion 
for  water  is  negative  thus  resulting  in  a  pressure  pulse,  which  has  a  rarefaction  preceding 
the  compressional  wave.  Two  other  observations  from  this  paper  are  worth  noting  for 
this  thesis.  One  is  the  formula  for  the  pressure  amplitude  expected  from  the 
instantaneous,  uniform  heating  of  a  long,  thin  rod-shaped  volume  of  water.  The 
expression  (attributed  to  Dolgoshein  and  Askarian)  is: 


p(r,e) 


K  E  c2 


sin 


71 L 
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sin# 


7lL  . 
—  sin# 
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\ 
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Eqn  1 


where  K  is  the  volume  coefficient  of  expansivity,  C  is  the  heat  capacity  of  water,  E  is 


the  total  energy  deposited,  r  is  the  distance  from  the  rod  (in  the  far  field),  c  is  the  velocity 
of  sound,  d  is  the  diameter  of  the  rod,  L  is  the  length  of  the  beam,  and  #  is  the  angle 
measured  relative  to  the  normal  to  the  rod.  Although  this  formula  is  specifically  for  the 
uniform  heating  of  a  rod,  a  maximum  will  occur  at  #  =  0°  for  any  linear  energy 
deposition  where  the  timescale  of  the  energy  deposition  is  fast  compared  to  the  acoustic 
timescale.  This  is  due  to  the  fact  that  the  contribution  to  the  pressure  from  all  points 
along  the  line  will  be  in  phase  along  the  acoustic  axis.  Therefore,  if  a  laser  pulse  were  to 
penetrate  deeply  into  the  water  creating  a  long  line  of  heating  with  an  exponential  decay, 
one  would  still  expect  to  get  a  pressure  maximum  nonnal  to  the  laser  beam.  Sulak  [1979] 
also  provides  an  equation  he  attributes  to  Bowen  for  the  pressure  amplitude  as  a  function 
of  time.  It  is: 


P(r,t) 


K  d  [dvw^J~r/c) 

4 nCp  dt  r 


Eqn  2 
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where  w(rj  )  is  the  heat  power  per  unit  volume  deposited  at  r  at  time  t,  dV  is  an 

infinitesimal  unit  of  the  heated  volume,  and  r  is  the  magnitude  of  the  distance  to  the  point 
where  the  pressure  is  measured.  From  this  formula,  it  can  be  seen  that  for  a  finite  laser 
pulse  in  the  time  domain,  the  pressure  depends  on  how  quickly  the  heat  power  is 
changing. 

The  other  observation  worth  noting  is  that  the  pulse  expected  from  this  heating  is 
a  bipolar  pulse  with  a  period  equal  to  the  time  required  for  the  pressure  wave  to  travel 
from  the  back  of  the  heated  volume  to  the  front  along  the  line  of  observation  plus  the 
energy  deposition  time. 

Maccabee  [1987]  also  derives  an  expression  for  the  pressure  expected  from 
heating  water.  His  treatment  is  different  from  the  one  we  see  in  Sulak,  because  he  looks 
specifically  at  a  situation  where  there  is  an  exponential  decay  in  intensity  along  the  laser 
beam.  The  intensity  of  the  laser  beam  in  water  at  depth  z  is  given  by  the  equation 

I(z )  =  I0e~az  Eqn  3 

where  a  is  the  attenuation  coefficient  of  water  and  70  is  the  laser  intensity  at  the  surface. 

Differentiating  both  sides  of  this  equation  with  respect  to  depth  yields 

=  -al0e~az  =  -«/(z) .  Eqn  4 

Since  the  change  in  intensity  of  the  beam  is  due  to  the  power  absorbed  by  the 
water,  the  power  absorbed  can  be  expressed  as: 

dE  dl  dl  T  (  \  ^  rz 

—  = - S  dz  = - dV  =  al(z)  AV  , 

dt  dz  dz 

where  S  is  the  cross-sectional  area  of  the  beam  and  AT  is  the  small  volume  in  which  the 
power  of  the  beam  is  deposited. 

This  absorbed  power  causes  both  the  temperature  and  volume  of  the  water  to 
increase  according  to 
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Eqn  5 


1  d2V  d2T 

AV  dt2  dt2  ’ 


where  p  is  the  coefficient  of  volume  expansion  and 


d2v 

dt 2 


is  the  acceleration  of  the 


volume  increase.  Furthermore,  the  change  in  temperature  with  respect  to  time  is  related 
to  the  heat  capacity  of  the  water,  the  volume,  and  the  power  absorbed  as: 


d2T  _  d 
dt 2  dt 


f  dE/  ^ 

/  dt 

PAVC„ 


Eqn  6 


where  p  is  the  density  of  water  and  Cp  is  the  specific  heat  per  unit  mass.  (It  is  assumed 

here  that  the  heat  capacity  and  density  of  the  water  do  not  change  during  the  heating 
process.)  Substituting  the  expression  for  power  above  yields: 


d  T  _  d 
dt2  dt 
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Eqn  7 


Rearranging  the  above  equations  yields  the  acceleration  of  the  volume  expansion 

al{z ) 


d  V  ad 
— —  =  aP  — 
dt  dt 


pcp 


AV  . 


Eqn  8 


This  is  not  the  expression  that  Maccabee  gives  in  his  paper.  He  omits  the  time 
derivative  and  gives  the  acceleration  of  volume  expansion  as: 


d2V  al(z ) 

■  =  ap — —  AF  . 


dr 


pcp 


Eqn  9 


This  is  clearly  an  error,  since  it  is  not  dimensionally  correct. 


Maccabee  goes  on  to  give  the  infinitesimal  contribution  to  the  pressure  from  the 
acceleration  of  the  volume  expansion  in  the  far-field  as: 

aBe~a' 

dp{t)  =  — — I0{t-r/c)dV,  Eqn  10 

4nrCp 

where  is  the  standard  retarded  time,  i.e.,  the  time  at  which  the  pressure 

contribution  from  dV  arrives  at  range  r .  Integrating  this  expression  over  the  heated 
volume  would  yield: 

p(r, t )  =  — — —  f al0 ( t-r  / c)e~a'Sdz  ,  Eqn  1 1 

\7irC  j 
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where  L  is  the  length  over  which  the  beam  deposits  its  energy.  This  expression  would 
be  identical  to  the  one  provided  in  Sulak’s  paper  if  the  derivative  with  respect  to  time  had 
not  been  dropped  from  the  formula  and  if  the  heat  power  deposited  per  unit  volume  had 

an  exponential  decay  with  distance,  i.e.,  w(t  -  =  al(t  -  y^j  =  ae~azI0(t  - 

Therefore,  the  correct  equation  for  Maccabee’s  case  should  be  given  as: 

p(r, t)  =  f i (t-r I c)e~azSdz .  Eqn  12 

4 nrC p  otJL 

This  is  an  important  correction,  because  it  predicts  that  the  pressure  produced 
depends  on  how  quickly  the  intensity  changes  with  time.  Assuming  that  the  heat 
capacity,  linear  expansion  coefficient,  and  attenuation  remain  constant,  an  increase  in 
pressure  can  be  attained  in  Type  1  sound  generation  either  by  increasing  the  energy  of  the 
laser  pulse  while  keeping  the  pulse  length  constant,  or  by  maintaining  the  pulse  energy 
while  decreasing  the  pulse  length. 

D.  TYPE  3  SOUND  PRODUCTION  -  SHOCK  WAVES 

Since  the  most  efficient  conversion  of  laser  power  to  acoustic  power  occurs  when 
the  energy  density  is  high  enough  to  produce  a  plasma,  this  type  of  sound  production  is 
more  interesting  in  terms  of  its  potential  as  an  acoustic  source.  Unfortunately,  under 
these  conditions,  the  resulting  acoustic  wave  falls  into  the  shock  regime.  Therefore,  the 
theory  describing  it  is  more  complicated.  Moreover,  the  factors  affecting  the  fraction  of 
laser  power  which  is  converted  to  an  acoustic  pulse  are  only  partly  understood.  The  best 
attempt  to  describe  the  complex  shock  wave  produced  by  a  high  energy  laser  pulse  is  that 
of  Vogel  et  al.  [1996,  1999].  The  motivation  of  their  work  is  to  understand  the  effects  of 
intraocular  laser  surgery  and,  in  particular,  to  find  a  regime  in  which  the  collateral 
damage  caused  by  the  expanding  bubble  and  shock  wave  is  minimized.  The  pulse 
energies  are  fairly  small;  however,  since  the  energy  densities  were  high  enough  to 
produce  Type  3  behavior,  their  treatment  is  relevant  to  this  thesis.  Of  course,  the  goal  of 
laser  acoustics  is  exactly  the  opposite  of  intraocular  surgery.  To  create  the  maximum 
possible  acoustic  signal,  laser  pulse  parameters,  which  maximize  the  transformation  of 
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laser  energy  into  mechanical  energy  are  desirable.  The  insights  gained  from  Vogel’s 
work  are  a  valuable  starting  point  for  laser  acoustics,  but  further  work  is  needed  to  ensure 
that  they  are  equally  valid  for  higher  energy  pulses. 

To  summarize  Vogel’s  insights  into  the  interaction  between  laser  pulses  and 
water,  a  timeline  of  the  interaction  will  be  presented  along  with  an  analysis  of  the  factors 
which  affect  the  efficiency  of  the  energy  conversion  from  laser  to  mechanical  energy. 
Although  they  investigated  the  process  for  several  different  wavelengths  ranging  from 
532  mn  to  1064  mn,  they  found  virtually  no  difference  in  the  behavior  as  a  function  of 
incident  wavelength.  This  is  not  too  surprising  since  the  laser  pulses  only  penetrated  the 
fluid  for  a  short  distance  before  reaching  the  focal  point.  Thus,  frequency  dependent 
absorption  was  not  an  issue.  This  would  certainly  not  be  the  case  for  an  application  in 
which  the  laser  pulse  penetrated  significant  distances  into  the  fluid  before  focusing. 

When  a  laser  pulse  enters  water,  the  absorption  coefficient  is  fairly  low  until  a 
threshold  energy  density  is  met,  which  causes  the  water  to  ionize  into  a  plasma.  Once  the 
plasma  is  formed,  the  density  of  free  electrons  keeps  the  light  absorption  extremely  high 
as  long  as  the  plasma  frequency  is  not  exceeded.  If  this  were  to  happen,  the  plasma 
would  become  perfectly  reflecting.  This  explains  the  substantially  higher  efficiency 
which  Maccabee  noted  for  the  Type  3  process.  The  intensity  first  becomes  high  enough 
to  form  a  plasma  at  the  focal  point.  Later,  as  the  pulse  intensity  rises,  the  plasma  grows 
toward  the  laser  beam  where  the  beam  cross-sectional  area  is  larger.  Since  the  cross- 
sectional  area  of  the  beam  increases  faster  for  larger  focusing  angles,  the  threshold 
intensity  required  to  ionize  the  water  can  only  be  achieved  close  to  the  focal  point.  For 
smaller  focusing  angles,  the  threshold  intensity  is  met  at  greater  distances.  This  results  in 
a  larger  initial  plasma  spot  size  and  lower  energy  density  than  the  larger  focusing  angle 
achieves.  This  concept  is  illustrated  in  Figure  2.  Here,  the  laser  beam  enters  from  the 
right.  The  apex  of  the  cone  is  the  focal  point  where  the  threshold  is  first  exceeded  before 
the  laser  pulse  reaches  its  maximum.  As  the  laser  pulse  builds  in  amplitude,  the  threshold 
is  reached  at  larger  and  larger  cross-sectional  areas  until  the  maximum  amplitude  of  the 
pulse  occurs.  For  the  larger  focusing  angle  on  the  right,  this  cross-sectional  area  is 
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achieved  closer  to  the  apex.  Since  the  volume  of  a  cone  is  given  by  —7ir"h ,  where  r  is 

the  radius  of  the  base,  the  plasma  volume  for  the  larger  focusing  angle  is  smaller.  The 
fact  that  the  plasma  is  growing  in  volume  towards  the  beam  tends  to  keep  electron 
densities  below  the  plasma  frequency  thus  contributing  to  the  ability  of  the  water  to 
absorb  the  incident  light  energy. 


volume 

Figure  2  Diagram  of  plasma  volume  for  small  (left)  and  large  (right)  focusing  angles. 

The  size  of  the  initial  plasma  spot  is  not  only  dependent  on  the  focusing  angle.  It 
also  depends  on  the  laser  pulse  length.  For  pulses  containing  the  same  total  energy  and 
shape,  the  peak  intensity  of  a  shorter  pulse  must  be  greater  than  the  peak  intensity  of  a 
longer  pulse.  This  results  in  a  larger  plasma  volume  since,  again,  the  intensity  reaches 
the  threshold  value  further  away  from  the  focal  point.  There  is  a  limit  to  the  extent  to 
which  the  pulse  can  be  lengthened  to  obtain  a  smaller  volume,  because  the  intensity  has 
to  exceed  the  threshold  in  order  to  get  sufficient  absorption.  Vogel  et  al.  [1999]  provide  a 
graph  of  the  percent  of  laser  energy  transmitted  through  the  fluid  as  a  function  of  the  total 
pulse  energy  for  a  30ps  pulse  length  and  various  focusing  angles.  This  graph  is  shown 
below.  At  lower  pulse  energies  transmission  is  very  high  since  the  threshold  has  not  been 
reached.  As  the  energy  exceeds  the  threshold  required  for  ionization,  the  transmission 
becomes  smaller  as  larger  amounts  of  the  laser  pulse  energy  are  absorbed  in  the  plasma. 
It  is  worthwhile  to  note  that  the  improvement  in  the  percent  of  energy  absorbed  is 
minimal  after  the  threshold  energy  has  been  exceeded  by  a  factor  of  about  50.  For  laser 
pulse  lengths  on  the  order  of  nanoseconds  and  below,  the  absorption  timeframe  is  fast  in 
comparison  to  the  rate  of  heat  conduction  and  shock  wave  propagation.  Therefore, 
immediately  following  the  laser  pulse,  the  laser  energy  absorbed  by  the  water  is  still 
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confined  to  a  volume  roughly  equal  to  the  initial  spot  size.  Since  the  laser  pulse 
timeframe  is  small  compared  to  the  consequent  thermodynamic  and  mechanical 
processes,  the  subsequent  evolution  is  relatively  insensitive  to  the  initial  laser  pulse 
length. 


Figure  3  Percent  transmission  versus  normalized  energy  for  a  30ps  pulse  at  various 
focusing  angles.  (From:  Vogel  et  al.  [1999]) 

Vogel  calculates  the  energy  required  to  heat  the  water  in  the  plasma  volume  to  the  boiling 
point  and  then  vaporize  it  as: 

Ev=PoVP[c(T2-Ti)  +  r]’  E9n13 

where  V  is  the  plasma  volume,  po  is  the  density  of  water,  c  is  the  heat  capacity  of 

water,  T. \  is  the  boiling  point,  7|  is  the  initial  temperature  of  the  water,  and  r  is  the 

specific  heat  of  vaporization  of  water.  No  mention  is  made  of  the  energy  required  to 
create  the  plasma  or  raise  it  to  its  final  temperature  at  the  end  of  the  laser  pulse. 
Presumably,  this  is  due  to  the  assumption  that  this  energy  is  completely  converted  to 
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mechanical  energy  in  the  form  of  the  bubble  expansion  and  shock  wave.  They  do 
investigate  the  energy,  which  is  lost  to  plasma  radiation  based  on  the  Stefan-Boltzmann 
Law.  However,  despite  the  high  plasma  temperature,  the  energy  lost  to  radiation  was 
found  to  be  negligible  in  all  cases. 

The  energy  density  within  the  plasma  volume  causes  it  to  heat  up,  vaporize,  and 
expand.  The  free  electrons  collide  and  transfer  their  energy  to  the  heavier  ions.  High 
speed  photographs  of  the  process  reveal  that  the  shock  front  detaches  from  the  plasma 
immediately  after  its  formation.  They  also  show  a  roughly  spherical  shock  front  and 
cavitation  bubble.  The  photograph  below  shows  a  long  cavitation  bubble  and  the  shock 
front  surrounding  it  produced  by  a  lmJ  30ps  pulse  in  the  upper  right.  In  comparison,  the 
lmJ  6ns  pulse  on  the  bottom  left  shows  a  shorter,  more  symmetrical  cavitation  bubble 
and  shock  front.  Assuming  that  the  size  of  the  cavitation  bubble  and  shock  front  are 
proportional  to  the  initial  spot  size,  this  photograph  thus  demonstrates  the  inverse 
relationship  between  plasma  spot  size  and  pulse  length. 


Figure  4  Photograph  of  cavitation  bubble  and  shock  wave  44ns  after  pulse.  (From: 
Vogel  et  al.  [1996]) 
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The  percentage  of  the  laser  pulse  energy,  which  is  converted  into  the  mechanical  energy 
of  the  cavitation  bubble  and  shock  wave,  depends  strongly  on  the  initial  energy  density. 
On  the  other  hand,  Vogel  found  that  the  ratio  of  the  energy  between  the  shock  wave  and 
the  bubble  were  remarkably  constant  at  about  2:1  over  a  wide  range  of  pulse  energies  and 
lengths  (from  a  50  pj  pulse  with  a  30ps  duration  to  lOmJ  and  6ns.)  Since  the  bubble 
energy  is  fairly  easy  to  measure  from  its  maximum  expansion  radius,  they  recommend  its 
use  in  determining  the  fraction  of  the  laser  pulse  energy  which  is  transformed  into 
mechanical  energy.  The  expression  used  for  the  bubble  energy  is  given  as  [Vogel  et  al. 
1999,  p.  273]: 

eb=-xKApo~  pv)>  Eqn14 

where  Rmax  is  the  maximum  radius  obtained,  po  is  the  hydrostatic  pressure  of  the  water 


the  bubble  expands  against,  and  pv  is  the  vapor  pressure  of  the  water.  They  use  a 

hydrostatic  pressure  of  O.IMPa  (1  atm)  corresponding  to  a  bubble  near  the  surface  of  the 
water  and  a  vapor  pressure  of  2330Pa  at  20°C.  Vogel  does  not  derive  this  formula,  but  it 
can  be  derived  from: 


W  =  fFfpdV. 


Eqn  15 


This  is  the  standard  formula  for  the  work,  W,  done  by  a  gas  expanding  against  a  constant 
external  pressure.  Direct  application  of  the  formula  would  yield: 

>Vr 


W=\r’pdV  =  po 


4  3 

-kRL-K 


V' 


Eqn  16 


J 


To  obtain  the  formula  they  give,  pVt  must  be  equated  with  pvVf  .  This  is  equivalent  to 


starting  with  an  ideal  gas  with  initial  volume,  Vj ,  at  hydrostatic  pressure  and  allowing  it 

to  expand  at  a  constant  temperature  until  the  external  pressure  is  equal  to  the  vapor 
pressure  of  water.  The  expansion  of  the  plasma  to  create  a  cavitation  bubble  is  hardly  an 
equilibrium  process,  nor  would  it  be  isothermal.  Moreover,  the  equation  of  state  for  the 
highly  condensed,  hot  plasma  probably  does  not  obey  the  ideal  gas  law.  However,  the 
correction  for  the  initial  volume  is  likely  to  be  small  for  cases  where  the  conversion  of 
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energy  to  mechanical  is  efficient.  At  the  very  least,  one  would  expect  the  maximum 
bubble  volume  to  be  proportional  to  the  energy  which  drives  the  bubble  expansion. 
Vogel  et  al.  [1999]  provide  a  graph  of  the  percent  of  the  laser  pulse  energy,  which  is 
converted  to  bubble  energy.  This  graph,  shown  below,  shows  the  dramatic  effect  of  laser 
pulse  length.  As  the  energy  is  raised  further  and  further  above  threshold,  there  is  little 
change  in  the  conversion  efficiency  at  any  of  the  pulse  lengths. 


0  50  100  150  200 

E/E* 

Figure  5  Percent  of  laser  pulse  energy  converted  into  cavitation  bubble  as  a  function 
of  pulse  length  and  nonnalized  pulse  energy.  (From:  Vogel  et  al.  [1999]) 


Due  to  the  difficulties  inherent  in  measuring  the  amplitude  of  the  shock  wave 
close  to  the  laser  pulse,  Vogel  and  his  collaborators  rely  on  an  interesting  interplay  of 
theory  and  experiment  to  determine  the  portion  of  energy,  which  goes  into  the  shock 
wave.  This  was  used  to  verify  the  2:1  proportioning  of  energy  between  the  shock  wave 
and  the  cavitation  bubble.  A  graph  of  the  shock  wave  as  a  function  of  time  is  shown 
below.  It  shows  the  classic  initial  shock  and  exponential  tail  of  the  shock  wave.  The 
width  of  the  shock  wave  appears  to  be  correlated  with  the  pulse  energy. 
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Figure  6  The  shock  wave  measured  10mm  from  the  plasma  volume  for  (a)  30ps, 
50  pj;  (b)  30ps,  lmJ;  (c)  6ns,  lmJ;  (d)  6ns,  lOmJ.  (From:  Vogel  et  al. 
[1996]) 


The  expression  used  by  Vogel  [1999]  for  the  energy  in  the  shock  wave  is: 


Es  = 


Anr~ 

A)co 


|  p2dt , 


Eqn  17 


where  p  is  the  nominal  density  of  water  and  CQ  is  the  nominal  sound  speed.  The 


pressure  can  be  measured  directly  by  a  small  hydrophone  with  a  quick  response  time  at 
distances  starting  at  about  lOnnn  from  the  laser  pulse.  At  shorter  ranges,  direct 
measurement  of  the  pressure  is  very  difficult.  Therefore,  they  used  a  couple  of  different 
techniques  to  try  to  estimate  the  pressure  closer  in.  The  first  technique  uses  the 
maximum  radius  of  the  cavitation  bubble,  the  laser  pulse  length,  and  the  plasma  volume 
to  predict  the  evolution  of  the  pressure  wave  in  the  surrounding  liquid.  Based  on  the 
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Gilmore  model  of  cavitation  bubble  dynamics,  calculations  were  performed  to  predict  the 
pressure  at  a  range  equal  to  six  times  the  initial  plasma  radius.  Using  the  pressure  from 
this  calculation  in  combination  with  the  equation  for  the  energy  in  a  shock  wave  above 

yields  the  shock  wave  energy,  Eflmo'e .  They  also  computed  the  shock  wave  energy 

using  a  different  technique  for  comparison  purposes.  The  other  technique  uses  the  energy 
left  10mm  from  the  center  (where  the  curvature  of  the  wavefront  is  large  enough  to  be 
measured  with  a  small  hydrophone)  and  the  energy  dissipated  as  the  shock  wave  travels 
from  its  initial  position  at  the  edge  of  the  plasma  volume  out  to  a  range  of  0.3mm.  The 
energy  remaining  at  10mm  is  referred  to  as  Esn 0mm,  and  the  energy  dissipated  up  to 


0.3mm  is  referred  to  as  EDiss .  The  equation  used  for  the  dissipation  of  energy  is  based 
on  the  Rankine-Hugoniot  equations.  The  increase  in  internal  energy  of  the  medium 
through  which  a  shock  wave  passes  is  given  by: 


1 


1 


3 


Ps(r)’ 


Eqn  18 


yPo  Ps{r)j 

where  p  is  the  density  of  the  shock  wave  and  p  is  the  pressure  of  the  shock  wave. 


Using  conservation  of  mass  at  the  shock  front  and  an  equation  of  state  they  attribute  to 
Rice  and  Walsh,  the  pressure  of  the  shock  wave  can  be  derived  from  the  velocity  of  the 
shock  front  as: 

Ps  =  ciPous  (10("'-C“)/C2  - 1)  +  P0  Eqn  19 

where  c0  is  the  normal  sound  of  sound  in  water,  ci  =  5190m/s,  C2  =  25,306m/s,  po  = 
998kg/m  ,  and  po  is  hydrostatic  pressure.  They  were  able  to  measure  the  velocity  of  the 


shock  front  photographically;  however,  as  the  velocity  approaches  the  nominal  speed  of 
sound,  the  uncertainty  of  the  pressure  calculation  becomes  large.  Therefore,  they  only 
used  this  method  to  calculate  the  dissipation  out  to  0.3mm.  From  the  increase  in  internal 
energy,  they  calculate  the  energy  dissipated  as: 
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Eqn  20 


n 

E Diss  =  j ^nr2ps  (r)As(r)  dr , 

'b 

where  ro  is  the  radius  at  the  onset  of  the  shock  wave  fonnation  and  l\  is  the  final  radius 

included.  The  density  at  the  shock  front  is  calculated  from  the  expression  for  the 
conservation  of  mass: 

UsPo={Us-Up)ps  Eqn 21 

and  the  expression  for  the  conservation  of  momentum: 

Ps-Po=usuPPo .  Eqn22 

where  u  is  the  particle  velocity  behind  the  shock  front.  They  claim  that  this  method 

underestimates  the  energy  in  the  shock  wave  since  it  fails  to  include  dissipation  from 
0.3mm  to  10mm.  However,  they  still  find  that  it  predicts  a  higher  shock  wave  energy  for 
all  cases  than  the  Gilmore  method.  They  attribute  this  discrepancy  to  the  possibility  that 
more  energy  is  lost  in  the  early  phases  of  the  shock  wave  -  before  the  Gilmore  method 
can  be  accurately  applied  -  than  in  the  0.3  -  lOnnn  region  that  is  omitted  in  the 
dissipation  method.  More  will  be  said  in  the  theory  chapter  about  potential  problems 
with  these  attempts  to  pin  down  the  shock  wave  energy. 

Other  conclusions  from  Vogel’s  work  of  interest  to  the  present  study  are  as 
follows: 

•  As  much  as  96%  of  the  incident  laser  energy  is  absorbed  by  the  plasma 
volume  for  a  6ns  lOmJ  pulse  [Vogel  et  al.  1999,  p.  275]. 

•  As  much  as  69%  of  the  absorbed  energy  goes  into  the  shock  wave  for  a 
6ns  lOmJ  pulse  [Vogel  et  al.  1999,  p.  275]. 

•  85-89%  of  the  shock  wave  energy  is  dissipated  in  the  first  0.2  -  0.3mm  of 
shock  wave  propagation  [Vogel  et  al.  1999,  p.  275]. 

•  Pressure  decays  as  r  2  until  the  pressure  falls  below  lOOMPa  [Vogel  et 
al.  1996,  p.  148]. 

•  After  the  pressure  falls  below  lOOMPa,  the  pressure  falls  as  approximately 
A1  [Vogel  et  al.  1996,  p.  148]. 
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To  the  extent  that  these  estimates  are  correct,  the  maximum  expected  efficiency  of 
the  electro-optic  to  acoustic  conversion  is  approximately: 

96%  x  69%  x  1 1  - 1 5%  =  7  - 1 0%  . 

This  represents  the  conversion  of  the  laser  pulse  into  useful  acoustic  energy  after  the 
strongly  dissipative  initial  stages  of  the  shock  wave  have  passed. 

Since  the  interest  in  this  study  is  on  the  pressure  amplitudes  obtainable  at  longer 
ranges,  it  is  also  useful  to  examine  more  closely  the  measured  shock  wave  parameters  at 
10mm.  Again,  using  the  later  (1999)  paper  and  the  more  efficient  6ns  pulse,  the  results 
can  be  summarized  in  the  following  table. 

Table  1  Experimental  data  measured  at  10mm  for  lmJ  and  lOmJ  laser  pulses  with 

length  of  6ns.  (From:  Vogel  et  al.  [1999])  Efficiency  computed  from  Vogel 
data. 


Laser  Energy 

lmJ 

lOmJ 

Pressure  (MPa) 

0.99 

2.62 

Duration  (ns) 

77 

148 

Energy  (pj) 

46.2 

622 

Efficiency  of  conversion 

5% 

6% 

(laser/acoustic  wave  at  10mm) 

Note  that  the  pressures  at  lOnnn  are  below  lOOMPa.  Therefore,  the  pressure  is 
expected  to  fall  as  J/  as  the  range  increases. 

E.  IDEAS  FOR  SHAPING  THE  ACOUSTIC  PULSE 

In  addition  to  Maccabee’s  early  suggestion  to  focus  the  laser  energy  on  an 
underwater  object,  there  are  two  other  innovative  ideas  increasing  the  efficiency  and/or 
directionality  of  the  laser  acoustic  source,  which  should  be  mentioned.  The  first  of  these 
techniques,  specifically  developed  with  airborne  laser  sources  in  mind,  seeks  to  maximize 
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the  acoustic  signal  through  the  placement  of  successive  pulses.  The  second  technique, 
recognizing  that  efficiency  is  proportional  to  the  energy  density  of  the  plasma,  seeks  to 
reduce  the  plasma  volume  through  non-linear  self-focusing  and  the  use  of  a  reverse  chirp 
signal. 

The  idea  of  moving  the  source  to  take  advantage  of  the  summation  of  successive 
pulses  was  studied  both  theoretically  and  experimentally  by  Bertholet  and  Busch-Visniac 
[1985].  They  recognized  that  by  moving  the  laser  spot  location  at  speeds  exceeding 
Mach  1,  the  wavefronts  from  successive  pulses  would  reinforce  each  other  and  create  a 
higher  amplitude  signal  than  could  be  achieved  by  a  single  pulse.  To  test  this,  they  used 
a  rotating  mirror  to  move  the  beam  across  the  surface  of  the  water.  They  also  modulated 
the  signal  by  using  a  Pockels  cell.  The  experimental  results  showed  a  maximum  increase 
of  20dB  in  the  pressure  amplitude  over  what  could  be  achieved  with  a  stationary  source. 
They  also  developed  an  ingenious  makeshift  modulation  technique  when  the  Pockels  cell 
failed.  By  laying  a  pattern  of  wooden  sticks  across  the  surface  of  the  water,  they  were 
able  to  modulate  the  laser  amplitude  as  it  scanned  the  water  surface  at  the  appropriate 
Mach  number  to  achieve  the  desired  wavefront  summation. 

More  recently,  Jones  et  al.  [2006]  reported  two  techniques,  which  could  be  used 
to  increase  the  electro-optic  to  acoustic  conversion  efficiency.  Noting  that  the  efficiency 
is  largely  dependent  on  the  energy  density  of  the  plasma  volume,  Jones’s  group  proposed 
using  two  different  techniques,  the  reverse  chirp  and  nonlinear  self-focusing  to  maximize 
the  energy  density  of  the  resultant  plasma. 

These  techniques  work  in  different  dimensions  to  decrease  the  size  of  the  initial 
laser  spot.  Nonlinear  self-focusing  (NSF)  relies  on  the  non-linearity  of  the  refractive 
index  at  high  intensity  to  provide  a  focusing  effect  in  the  water  medium.  This  results  in  a 
constriction  in  the  transverse  dimension.  At  the  same  time,  a  constriction  along  the 
longitudinal  dimension  is  achieved  by  taking  advantage  of  the  dispersive  properties  of 
water.  By  starting  out  with  a  negatively  chirped  laser  pulse,  the  group  velocity  dispersion 
of  water  causes  a  longitudinal  compression  as  shown  in  the  figure  below.  Unfortunately, 
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the  beam  has  to  travel  a  relatively  large  distances  to  achieve  the  required  focusing,  i.e., 
about  8m  for  the  experiments  conducted  by  Jones.  This  results  in  some  attenuation  in  the 
pulse  energy  based  on  the  wavelength  of  the  laser  pulse  employed. 


Relatively  low  intensity 

High  intensity 

negatively  chirped  laser  pulse 

| 

compressed 

— — 

Longitudinal 

pulse 

◄ - ► 

pulse  duration:  T(z) 

compression 

Figure  7  Transverse  and  longitudinal  focusing  from  a  negatively  chirped  laser  pulse 
in  water  (From:  Jones  et  al.  2006) 

F.  EXPERIMENTAL  RESULTS 

Several  factors  that  affect  the  efficiency  of  the  electro-optic  to  acoustic 
conversion  were  discussed  in  the  previous  sections.  At  a  general  level,  efficiency  is 
affected  by  pulse  length,  focusing  angle,  total  pulse  energy,  and  wavelength  (for  pulses 
that  penetrate  into  the  water.)  To  get  a  quick  overview  of  the  pressure  amplitudes,  which 
were  achieved  by  various  research  groups,  the  peak  amplitudes  from  their  “best  shots” 
are  presented  in  the  graph  below  as  a  function  of  pulse  energy.  In  order  to  make  the 
comparisons  meaningful,  the  pressure  amplitude  is  extrapolated  to  a  range  of  lm.  This  is 
justified  in  that  all  measurements  were  taken  for  pressure  amplitudes  below  lOOMPa,  the 

transition  point  between  y2  and  /  spreading.  Following  the  graph,  a  table  is 
provided  to  detail  the  parameters,  which  resulted  in  the  plotted  measurement. 
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Peak  pressure  (Pa) 


Peak  pressure  vs  laser  pulse  energy 


fi 


Pulse  energy  (J) 

Figure  8  Peak  pressures  achieved  by  various  research  groups.  Values  are 
extrapolated  to  a  range  of  lm  using  J/  spreading. 
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Table  2  Parameters  used  to  achieve  peak  pressures  shown  in  Figure  8 


Group 

Pulse 

energy 

Pulse 

length 

Spot  Size 

Focusing 

Comments 

Bertholet  and 

Busch- 

Visniac 

[1985] 

1.5J 

1 .2ms 

0.8cm2 

not  focused 

1.38Pa  at  0.5m 

1 .06  p  m  wavelength 

a  =  13.7m1 

no  shock  wave 

Maccabee 

[1987] 

1.6J 

200ns 

0.2cm2 

focused  on 
surface 

2.3MPa  at  4.9cm 
directly  below 

10.6  pm  wavelength 

“Type  2”  - 
evaporative 

Vogel  et  al. 

[1996] 

lOmJ 

6ns 

45  pm2 

focused  -  22° 
convergence 
angle 

2.62MPa  at  10mm 

1 .06  p  m  wavelength 

Jones  et  al. 

[2003] 

120J 

5  ns 

0.2  -  0.5cm2 

focal  length 
of  2m 

1 8MPa  at  1 3mm 

527nm  wavelength 

a  =  0.02m1 

focused  on  plastic 
foil 

Jones  et  al. 

[2006] 

2.8mJ 

50fs 

unknown 

unknown 

49kPa  at  5mm 

400nm  wavelength 

a  =  0.1m'1 

Since  the  peak  intensity  of  the  acoustic  signal  should  be  roughly  proportional  to 
the  energy  in  the  laser  pulse,  it  is  also  interesting  to  fonn  a  relative  figure  of  merit  (FOM) 
as  the  peak  pressure  squared  divided  by  the  pulse  energy.  This  FOM  fails  to  take  into 
consideration  the  duration  either  of  the  acoustic  pulse  or  any  asymmetries  in  the  radiation 
pattern.  However,  these  factors  were  not  provided  by  all  groups.  Therefore,  the  FOM 
defined  in  this  way  cannot  be  used  as  a  precise  measure  of  the  conversion  efficiency. 
Table  3  shows  the  FOM  calculated  from  the  data  above. 
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Table  3  Ratio  of  peak  pressure  squared  to  laser  pulse  energy  achieved  by  various 
research  groups. 


Group 

r  2  /  \ 

FOM  Ppk/E 

y  ^  laser  y 

Bertholet  and  Busch-Visniac 

[1985] 

0.32  PaYj 

Vogel  et  al. 

[1996] 

6.9x\06PaY 

Jones  et  al. 

[2006] 

2.\x\07  PaY 

Jones  et  al. 

[2003] 

4.8x10  *Pa2/j 

Maccabee 

[1987] 

7.6x109jP^ 

From  the  FOM  data  shown  above,  clearly  the  result  from  Bertholet  and  Busch- 
Visniac  shows  the  smallest  conversion  efficiency  from  electro-optic  to  acoustic  energy. 
This  is  consistent  with  a  Type  1  heating  process  resulting  from  the  long  pulse  length  and 
consequently  low  peak  laser  intensity.  Maccabee’s  1.6J  shot  shows  the  highest 
conversion  efficiency  even  though  he  considers  it  to  be  the  result  of  a  Type  2  evaporative 
process  on  the  surface  of  the  water.  This  is  a  surprising  result  and  raises  the  question  of 
whether  the  process  was  purely  Type  2,  or  whether  there  was  plasma  formation  as  well. 
The  2003  shot  of  Jones  et  al.  shows  that  a  much  higher  energy  pulse  (120J)  can  also 
achieve  very  high  conversion  efficiencies  when  the  plasma  spot  size  is  minimized 
through  the  use  of  a  foil  target.  This  group’s  2006  shot  shows  just  an  order  of  magnitude 
lower  conversion  efficiency  with  a  much  shorter  pulse  length  and  no  foil.  Finally,  Vogel 
et  al.  achieved  a  conversion  efficiency  about  an  order  of  magnitude  lower  than  Jones 
[2006],  This  is  somewhat  surprising  since  their  beam  was  focused,  and  they  used  a 
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longer  pulse  length.  According  to  Vogel,  both  of  these  factors  should  have  contributed  to 
the  smaller  spot  size  and  higher  initial  energy  density  associated  with  a  high  conversion 
rate. 
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III.  THEORETICAL  ACOUSTIC  PULSE  ESTIMATES 


The  physics  of  the  interaction  between  a  high  energy  laser  pulse  and  water  is  not 
understood  well  enough  to  calculate  the  exact  details  of  the  acoustic  wave  produced 
based  solely  on  the  laser  pulse  parameters,  placement  of  the  pulse,  and  physical 
properties  of  the  water.  It  is  not  even  clear  that  the  problem  can  be  well  bounded  due  to 
the  number  of  factors  involved.  However,  it  is  possible  to  develop  some  intuition  into  the 
factors,  which  affect  the  characteristics  of  the  acoustic  pulse.  Therefore,  the  first  section 
of  this  chapter  lists  all  the  factors  either  which  have  been  shown  or  which  are  considered 
likely  to  play  a  part  in  detennining  the  characteristics  of  the  resultant  acoustic  pulse. 
These  factors  are  then  discussed  in  an  attempt  to  clarify  the  roles  they  play  in  determining 
the  acoustic  pulse.  The  last  section  of  this  chapter  looks  at  several  of  Vogel’s 
calculations  in  more  detail  and  speculates  qualitatively  about  the  effect  of  the 
assumptions  made  on  the  conclusions  drawn  from  that  work. 

A.  FACTORS  AFFECTING  THE  ACOUSTIC  PULSE  PRODUCED  BY  A 

HIGH  ENERGY  LASER  PULSE 

In  general,  the  parameters  which  affect  the  acoustic  signal  produced  by  a  high 
energy  laser  pulse  can  be  divided  into  two  basic  categories,  i.e.,  the  laser  pulse 
characteristics  and  the  characteristics  of  the  acoustic  medium.  In  some  cases  the  effect  of 
these  characteristics  are  inextricably  linked. 

1.  Laser  Pulse  Characteristics 

The  laser  pulse  characteristics,  which  detennine  the  resultant  acoustic  signal,  can 
be  summarized  as: 

•  Total  energy  per  pulse 

•  Pulse  length  in  time  domain 

•  Pulse  shape  in  time  domain 

•  Pulse  cross-section 

•  Focus  angle 
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•  Wavelength 

•  Angle  of  incidence  into  medium 

•  Pulse  repetition  rate 

While  the  total  energy  of  the  laser  pulse  sets  an  upper  limit  to  the  energy  which 
ends  up  in  the  resultant  acoustic  signal,  other  factors  can  cause  the  percent  energy 
conversion  to  vary  from  almost  0%  to  a  high  of  about  10%.  For  example,  operating  at  a 
wavelength,  which  suffers  little  absorption  in  water,  causes  the  energy  transmission 
through  the  water  to  be  very  high.  Thus,  little  of  the  incoming  electro-optic  energy  is 
converted  into  an  acoustic  pulse.  On  the  other  hand,  if  the  energy  density  is  high  enough 
to  form  a  plasma,  more  of  the  electro-optic  energy  is  absorbed.  The  resulting  heating  and 
vaporization  of  the  water  transfers  energy  to  the  acoustic  wave.  This  plasma  formation  is 
the  key  to  obtaining  high  efficiencies  since  even  wavelengths  which  absorb  well  in  water 
will  not  deposit  sufficient  energy  in  a  small  enough  volume  to  attain  good  conversion 
efficiencies.  In  general,  the  greater  the  energy  density  in  the  water  immediately 
following  the  laser  pulse,  the  greater  the  efficiency  of  the  electro-optic  to  acoustic 
conversion.  For  a  given  total  pulse  energy,  the  peak  intensity  of  the  laser  pulse  is 
inversely  proportional  to  its  duration.  Therefore,  smaller  pulse  lengths  in  the  time 
domain  tend  to  result  in  higher  efficiencies.  On  the  other  hand,  if  the  energy  density  is 
too  high,  electron  densities  may  be  high  enough  to  exceed  the  plasma  frequency.  In  this 
case,  all  remaining  electro-optic  energy  would  be  reflected.  Thus  the  upper  limit  to  the 
energy  density  is  detennined  by  the  point  where  the  plasma  frequency  is  reached.  All 
other  factors  being  equal,  it  might  be  possible  to  harness  more  of  the  energy  of  a  pulse  by 
shaping  it  in  such  a  way  that  a  plasma  is  formed  very  early  so  that  absorption  is  high 
throughout  the  following  pulse  evolution.  The  changing  speed  of  light  in  the  water  as  a 
function  of  the  heating,  vaporization,  ionization,  and  pressure  is  another  issue,  which 
affects  the  shape  of  the  electro-optic  pulse  as  it  travels  though  the  water.  Jones  et  al. 
[2006]  have  suggested  how  the  natural  dispersion  of  the  medium  can  be  harnessed  to 
increase  the  energy  density  of  a  reverse  chirp  pulse.  In  general,  however,  refraction  of 
the  laser  light  away  from  regions  of  lower  density/higher  temperatures  and  scattering  off 
bubbles  has  not  been  considered.  It  is  likely  that  these  effects  would  be  more  important 
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for  multiple  pulses  directed  in  rapid  succession  at  the  same  area.  Since  multiple  pulses 
will  be  required  for  achieving  low  frequency  acoustic  signals,  refraction  and  scattering 
issues  may  be  important. 

The  cross-section  of  the  laser  beam  also  plays  an  important  role  in  the  energy 
density  attained.  The  cross-section  of  the  beam  in  the  water  depends  on  the  original 
cross-sectional  characteristics  of  the  beam  as  well  as  the  distance  to  the  water, 
atmospheric  conditions  affecting  the  beam  propagation,  difference  in  refractive  index  of 
water  and  air,  beam  focusing,  and,  again,  as  pointed  out  by  Jones  et  al.  [2006]  the  non¬ 
linear  self-focusing  effect.  Since  the  efficiency  of  the  conversion  is  dependent  on 
whether  the  energy  density  is  high  enough  to  create  a  plasma,  efficiency  is  favored  by 
having  a  cross-sectional  area  small  enough  to  attain  the  threshold  energy  density  required 
for  ionization.  External  focusing  with  large  focusing  angles  led  in  Vogel’s  work  to 
higher  efficiencies  since  the  cross-sectional  area  required  to  reach  the  threshold  energy 
density  occurred  closer  to  the  focal  point.  This  resulted  in  a  higher  overall  energy  density 
of  the  pulse. 

The  amount  of  the  original  pulse  energy,  which  enters  the  water,  depends  on  the 
transmission  coefficient  of  the  air/water  interface.  This  coefficient  depends  on  the  angle 
of  incidence  and  the  indices  of  refraction  for  water  and  air.  Of  course,  it  also  depends  on 
how  smooth  the  surface  is.  For  ocean  applications,  the  surface  of  the  ocean  might  be 
locally  smooth  in  the  vicinity  of  the  laser  beam,  but  the  angle  of  incidence  would  be 
fairly  random.  On  the  other  hand,  since  the  speed  of  light  is  slower  in  water  than  it  is  in 
air,  there  is  no  critical  angle  at  which  the  transmission  coefficient  goes  to  zero. 

2.  Characteristics  of  Acoustic  Medium 

The  characteristics  of  the  acoustic  medium,  in  this  case  water,  which  affect  the 
ability  of  the  medium  to  transform  the  electro-optic  energy  from  the  laser  pulse  into 
acoustic  energy  can  be  summarized  as: 

•  Speed  of  sound  (which  is  a  function  of  temperature,  pressure,  and  salinity) 

•  Density 

•  Viscosity 
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•  Heat  capacity 

•  Heat  conduction 

•  Heat  of  Vaporization 

•  Ionization  energy 

•  Shape  and  location  at  which  electro-optic  energy  is  deposited 

As  the  laser  energy  is  absorbed  by  the  water,  it  heats  up.  If  the  timeframe  of  the 
laser  pulse  is  fast  compared  to  the  time  it  takes  for  heat  to  conduct  away  or  for 
vaporization  to  occur,  the  spot  where  the  energy  is  deposited  can  be  envisioned  as 
superheated,  ionized  water.  The  temperature  of  the  water  is  determined  by  the  heat 
capacity  of  the  water  at  constant  volume  and  by  the  energy  lost  to  ionization  and  heat 
conduction  over  the  timeframe  of  the  laser  pulse.  This  superheated  water  is  going  to 
expand.  Vogel’s  experimental  evidence  shows  almost  no  loss  of  energy  in  the  form  of 
radiation.  This  is  probably  due  to  the  density  of  the  water  in  the  laser  spot.  Collisions 
between  electrons  and  ions  most  likely  lead  to  a  quick  recombination  and  the  consequent 
transfer  of  the  energy  to  heat.  This  fast  conversion  of  the  energy  lost  to  ionization  back 
to  heat  energy  would  also  explain  why  Vogel  did  not  include  any  mention  of  ionization 
energy  in  his  energy  balance. 

The  expanding  spot  creates  a  shock  wave,  which  detaches  almost  immediately 
from  the  bubble  of  vaporized  water.  The  bubble  will  continue  expanding  until  its  internal 
pressure  is  balanced  by  the  external  pressure.  It  will  then  oscillate  with  damping  and 
gradually  get  reabsorbed  in  the  water  as  the  temperature  decreases. 

The  details  of  the  shock  wave  propagation  depend  on  the  speed  of  sound  in  water 
at  the  various  local  temperatures  and  densities  as  well  as  the  small  viscosity  of  the  water. 
Energy  is  lost  to  the  shock  wave  due  to  heat  conduction.  The  placement  of  the  pulse 
relative  to  the  surface  of  the  water  plays  an  important  role  in  the  received  acoustic  signal 
since  it  reflects  the  acoustic  pulse  with  a  180°  phase  shift. 
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B.  ENERGY  BALANCE 

Since  it  is  impossible  to  predict  the  precise  acoustic  characteristics  expected  based 
solely  on  the  characteristics  of  the  input  laser  pulse  and  the  water,  an  understanding  of  the 
energy  balance  of  the  process  is  particularly  valuable.  In  this  regard,  several  of  the 
equations  used  by  Vogel  et  al.  [1996,  1999]  bear  further  examination. 

To  start  with,  Vogel  et  al.  [1999]  give  the  following  figures  for  the  energy  balance 
of  their  6ns  lOmJ  pulse: 

•  Energy  reflected  off  laser  spot  -  0.8% 

•  Energy  scattered  off  laser  spot  -  0.5% 

•  Energy  transmitted  though  laser  spot  -  2.6% 

•  Energy  absorbed  by  water  -  96. 1  % 

•  Mechanical  energy  needed  to  form  bubble  -  29.4% 

•  Shock  wave  energy  -  49.3  -  68.5%  (depending  on  how  measured) 

•  Energy  of  V aporization  -6.5% 

•  Energy  radiated  -  6  x  10 '4% 

•  Energy  unaccounted  for  -  1 .3% 

The  energy  reflected  off,  scattered  by,  and  transmitted  through  the  laser  spot  are 
measured  by  calibrating  a  receiver  by  assuming  that  100%  of  the  light  incident  on  the 
sample  at  low  energies  is  transmitted.  Thus,  these  measurements  do  not  include  energy 
lost  at  any  of  the  optical  interfaces  other  than  the  laser  spot  itself. 

The  energy  of  vaporization  is  calculated  using  the  equation: 

Ev  =  P0Vp  [c  (T2  -  7] )  +  r]  ,  Eqn  23 

where  poV  is  the  volume  of  the  laser  spot,  c  is  the  heat  capacity,  and  r  is  the  latent  Heat 

of  Vaporization.  As  noted  previously,  this  equation  does  not  account  for  the  energy 
required  to  attain  temperatures  above  the  boiling  point  of  water  at  atmospheric  pressure, 
100°C.  It  assumes  that  only  the  water  in  the  volume  of  the  initial  laser  spot  is  vaporized. 
They  use  the  Heat  of  Vaporization  and  the  heat  capacity  for  isobaric  processes.  The 
assumption  that  the  heating  and  vaporization  are  isobaric  processes  and  the  use  of  values 
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measured  for  equilibrium  processes  at  specific  temperatures  are  not  accurate  for  this  non¬ 
equilibrium  process.  Furthermore,  since  the  temperature  of  the  bubble  is  likely  to  be 
considerably  higher  than  100°C,  this  approach  is  expected  to  underestimate  the  energy 
required  to  vaporize  and  heat  the  water. 

The  energy  required  to  form  the  cavitation  bubble  is  calculated  using  the 
equation: 

EB  =  27rRL(Po~Pv)’  Eqn24 

where,  again,  Rnmx  is  the  maximum  radius  obtained,  pQ  is  the  hydrostatic  pressure  of  the 
water  the  bubble  expands  against,  and  pv  is  the  vapor  pressure  of  the  water.  Again,  the 
issue  with  the  use  of  this  formula  is  that  using  dW  =  pdV  to  calculate  the  work  done  by 
an  expanding  gas  assumes  an  equilibrium  process. 

The  energy  of  the  shock  wave  is  also  a  difficult  calculation  to  do  under  the 
conditions  of  a  non-equilibrium  process.  Vogel  does  this  calculation  in  two  different 
ways  -  both  of  which  admittedly  underestimate  the  energy  in  the  shock  wave.  One  uses 
the  measurement  of  the  pressure  at  a  distance  of  10mm  to  calculate  the  energy  left  at  this 
point.  That  quantity  is  added  to  the  energy  dissipated  up  to  a  range  of  0.3mm.  This 
method  is  considered  to  underestimate  the  initial  energy  in  the  shock  wave  because 
dissipation  from  0.3mm  to  10mm  is  not  included.  The  other  method  used  to  estimate  the 
energy  in  the  shock  wave  is  to  use  the  equation: 


Since  the  pressure  cannot  be  measured  directly  close  to  the  laser  spot,  the  pressure  was 
estimated  using  the  Gilmore  model  and  experimentally  determined  values  for  the 
maximum  bubble  radius,  laser  spot  size,  and  pulse  duration.  Aside  from  any  questions  on 
the  ability  of  the  Gilmore  model  to  predict  the  shock  wave  pressure,  it  should  be  pointed 
out  that  the  equation  used  above  for  the  energy  in  a  shock  wave  is  derived  from  linear 
acoustics.  It  comes  from  the  fact  that  the  instantaneous  intensity,  l(t),  of  an  acoustic 
wave  is  given  by: 
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I(t)  =  p(t)u(t), 


Eqn  26 


where  u(t )  is  the  particle  velocity.  For  a  linear,  spherical  wave,  the  magnitude  of  the 


pressure  and  particle  velocity  is  given  by  KFCS  by: 


and 


pit )  =  —cos  (cot-kr) 


l  A  l 

u(t)  = - cos  (cot-kr-  6*) , 


poc  r  cos  9 


Eqn  27 

Eqn  28 


where  —  is  the  amplitude  of  the  pressure  at  range  r,  co  is  the  angular  frequency,  k  is  the 
r 

wavenumber,  and  0  is  the  phase  difference  between  the  pressure  and  particle  velocity. 
For  a  spherical  wave,  this  phase  difference  depends  on  range  and  is  given  by: 

tan-1  —  .  Eqn  29 

kr 

This  is  significant,  because  at  very  small  values  of  kr,  the  pressure  and  particle  velocity 
are  close  to  90°  out  of  phase.  Thus,  high  particle  velocities  yield  little  in  the  way  of 
acoustic  pressure  for  small  acoustic  sources.  Integrating  the  instantaneous  intensity  over 
one  period  yields  the  expression: 

P2 

I  = - ,  Eqn  30 

2PoC 

where  P  is  the  pressure  amplitude  at  range  r.  For  a  finite  wave  containing  different 
frequencies,  each  frequency  component  would  behave  the  same,  thus  yielding  the 
expression  given  by  Vogel.  However,  the  relationship  between  pressure  and  particle 
velocity  is  incorrect  in  the  highly  non-linear  regime  of  shock  waves.  The  relationship 
between  pressure  and  particle  velocity  given  by  Vogel  as: 

Ps  =  ci Pa11.,  (lO("’-")/C2  -l)  +  P0  Eqn  3 1 

should  be  used  instead  to  estimate  the  energy  in  the  shock  wave. 
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IV.  EXPERIMENTS 


Experiments  were  conducted  with  Ted  Jones  at  the  Naval  Research  Laboratory, 
Code  6795  during  the  week  of  July  23-27,  2007,  using  a  Tiisapphire  laser.  The  goal  of 
these  experiments  was  to  get  a  preliminary  understanding  of  the  equipment  and 
procedures  involved  in  laser-generated  acoustics  as  well  as  to  do  some  rough 
measurements  of  the  acoustic  signal  generated  by  a  fairly  low  energy  pulse.  The 
experiments  were  successful  in  measuring  the  pressure  and  its  range  dependence  starting 
at  a  range  of  a  few  centimeters  from  the  source.  Experiments  were  also  conducted  to  get 
a  feel  for  the  directionality  of  the  acoustic  source;  however,  these  were  inconclusive  due 
to  the  possible  role  of  surface  interference  in  determining  the  received  pressure 
amplitudes. 

The  pulses  used  in  these  experiments  had  an  approximate  pulse  length  of  lOOfs,  a 
400nm  wavelength,  and  an  approximate  energy  of  2mJ.  They  were  directed  into  an  open 
air  water  tank  that  had  a  volume  of  lm  .  The  tank  was  approximately  half  full  and  the 
pulse  was  directed  at  the  surface  of  the  water  using  mirrors  and  lenses.  A  diagram  of  the 
optical  set-up  is  shown  in  Figure  9,  and  a  photograph  of  the  equipment  is  shown  in  Figure 
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Mirrors 


Laser  pulse 


Figure  9  Diagram  of  Experimental  Setup. 
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Figure  10  Picture  of  experimental  setup. 

A  negatively  chirped  laser  pulse  was  not  used  in  any  of  these  experiments  so  no 
longitudinal  compression  of  the  laser  beam  was  expected  once  the  laser  hit  the  water. 
Also,  because  the  tank  was  not  deep  enough  for  nonlinear  self  focusing  (NSF)  to  occur,  a 
lens  was  used  to  focus  the  laser.  The  lens  was  moved  up  and  down  in  an  attempt  to  focus 
the  laser  just  below  the  surface  of  the  water,  but  without  the  use  of  a  high  speed  camera  it 
was  very  difficult  to  tell  with  the  naked  eye  if  the  laser  spot  was  on  or  below  the  surface 
of  the  water.  Instead,  the  lens  was  moved  until  a  maximum  response  was  observed  on  the 
oscilloscope.  A  telescope  was  used  to  make  the  cross  section  of  the  beam  larger  to 
prevent  energy  loss  from  the  beam  and  damage  to  the  mirrors  and  lenses  due  to 
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filamentation.  Filamentation  is  caused  by  self  focusing  of  the  laser  beam  which  results  in 
intensities  high  enough  to  produce  a  plasma.  The  plasma  formation  in  air  removes 
energy  from  the  laser  beam  before  it  reaches  the  water. 

Acoustic  data  was  obtained  using  a  low  frequency  hydrophone  and  two  high 
frequency  hydrophones.  The  low  frequency  hydrophone  was  a  TC-4034  manufactured 
by  Reson.  The  hydrophone  was  connected  to  the  oscilloscope  through  a  VP  1000 
preamplifier  also  manufactured  by  Reson.  The  high  frequency  hydrophones  were  1mm 
needle  hydrophones  manufactured  by  Precision  Acoustics  LTD.  They  were  connected  to 
the  oscilloscope  through  HA60  preamplifiers  also  manufactured  by  Precision  Acoustics 
LTD.  Using  the  voltage  from  the  oscilloscope,  V,  the  pressure,  P,  was  determined  using 
the  sensitivity,  M,  of  the  hydrophones  from  the  hydrophones’  calibration  curves  provided 
in  the  Appendix  along  with  the  preamplifier  gain,  G,  in  dB. 

P  = - Eqn  32 

M(  1(P) 

A.  DIRECTIONALITY  STUDY 

The  first  experiment  conducted  was  an  attempt  to  determine  the  directionality  of 
the  acoustic  pulse  from  a  laser-generated  acoustic  source  using  the  high  frequency 
hydrophone.  The  received  pulse  is  shown  in  Figure  11.  It  can  be  seen  in  this  figure  that 
the  period  of  the  pulse  is  about  l.lps,  which  leads  to  a  pulse  frequency  of  about  0.91 
MHz. 
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Figure  1 1  Pulse  shape  from  high  frequency  hydrophone  achieved  during  the  30° 
directionality  trial. 

Since  there  was  no  high  speed  camera  available  to  photograph  the  interaction  of 
the  laser  pulse  with  the  water,  it  was  impossible  to  measure  the  precise  distance  between 
the  laser  spot  and  the  hydrophone.  Therefore,  the  position  of  the  hydrophone  was 
adjusted  until  the  elapsed  time  between  the  initiation  of  the  laser  pulse  and  the  arrival  of 
the  acoustic  pulse  at  the  hydrophone  was  approximately  the  same  at  each  angle.  The 
angle  itself  had  to  be  estimated  with  a  protractor  by  assuming  the  location  of  the  laser 
spot  was  on  the  surface  of  the  water  along  the  beam  direction.  Care  was  taken  to  point 
the  hydrophone  as  close  to  directly  at  the  laser  spot  as  possible  while  using  the  naked  eye. 
However,  based  on  the  relatively  flat,  1MHz  line  on  the  directionality  plot  for  the  high 
frequency  hydrophone  in  the  Appendix,  the  pointing  direction  of  the  hydrophone  should 
not  make  a  noticeable  difference.  A  basic  diagram  is  shown  in  Figure  12. 
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Figure  12  Diagram  of  basic  setup  for  directionality  experiment. 

The  results  are  shown  in  Figure  13  with  0°  being  directly  below  the  laser  spot. 
Each  angle  shows  the  data  from  three  laser  pulses  except  for  45°  where  the  data  from  two 
of  the  pulses  was  corrupted.  The  hydrophone  was  not  moved  between  pulses  at  each 
angle.  The  peak  pressures  attained  from  these  trials  were  higher  at  0°  and  30°  and  then 
fell  off  at  sharply  at  45°  and  above.  This  observation  might  be  consistent  with  a 
directional  pulse  if  the  beam  entered  the  water  at  an  angle  such  that  the  mainlobe  of  the 
beampattem  occurred  around  30°,  and  if  the  pulse  energy  were  deposited  along  a  line 
whose  length  was  larger  than  its  diameter.  However,  the  pulse  was  focused  onto  a  small 
spot,  and  the  direction  of  entry  was  nominally  normal  to  the  surface.  Therefore,  it  is 
more  likely  that  the  observed  drop  in  peak  pressure  is  due  to  interference  between  the 
direct  acoustic  pulse  and  its  reflection  off  the  surface. 
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Directionality  Study  Using  High  Frequency  Hydrophone 
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Figure  13  Directionality  study  of  laser-generated  acoustic  source  using  a  high 
frequency  hydrophone. 


The  effect  of  surface  interference  can  be  explained  by  an  analysis  of  the  diagram 
shown  below. 
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Figure  14  Surface  interference  geometry  for  calculating  the  path  length  difference 
between  direct  and  surface  reflected  path. 


The  laser  pulse  is  deposited  at  a  depth  of  d  beneath  the  surface  and  a  hydrophone 
is  located  at  a  range  r  and  an  angle  0  from  it.  The  surface  reflection  travels  a  longer 
distance,  rs ,  to  the  hydrophone.  Using  the  Pythagorean  Theorem,  the  length  of  the 

reflected  path  can  be  expressed  in  terms  of  the  source  depth,  hydrophone  separation,  and 
angle  as: 

r2  =r2sin2#  +  (2d  +  rcos#)2.  Eqn  33 

Multiplying  this  out  and  simplifying  leads  to: 

r2  =  r2  +  4rd cos 6  +  4d2 .  Eqn  34 

When  d  is  much  less  than  r  then  this  can  be  further  simplified  to: 


r~  =  r  +  4 rd  cos  6  -  r~ 


1  + 


4  d  cos  0 


Eqn  35 


Since  the  cosine  has  a  maximum  value  of  one  and  d  is  assumed  to  be  small,  the 
length  of  the  reflected  path  can  be  approximated  by: 

2d  cos  6 


r  =  r 


1  +  - 


Eqn  36 
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Thus,  the  path  length  difference  between  the  direct  and  reflected  path  can  be 
approximated  as: 

rs  -r  =  2d cos#  .  Eqn  37 

The  maximum  path  length  difference,  2d  , occurs  when  9=0°  which  is  when  the 
hydrophone  is  directly  below  the  laser  spot.  The  minimum  path  length  difference  occurs 
when  0  =  90°,  or  when  the  hydrophone  is  perpendicular  to  the  laser  direction.  Since  the 
direct  and  reflected  path  lengths  are  180°  out  of  phase  due  to  the  phase  shift  resulting 
from  the  air-water  impedance  difference,  the  two  paths  interfere  destructively  if  they 
arrive  simultaneously.  With  the  measured  temporal  pulse  length  of  1 . 1  p  s,  the  physical 
length  of  the  pulse  (assuming  a  propagation  speed  of  1480m/s)  is  about  1.6mm.  Path 
length  differences  less  than  about  a  millimeter  would  allow  the  negative  part  of  the 
reflected  pulse  to  arrive  during  the  compressional  part  of  the  direct  path  pulse.  For 
example,  if  the  laser  spot  were  focused  at  a  depth  of  1mm,  the  path  length  between  the 
direct  and  reflected  path  would  be  about  1mm  at  an  angle  of  60°.  At  smaller  angles,  the 
path  length  difference  would  be  larger,  and  the  pulses  would  start  to  separate.  Therefore, 
the  fact  that  the  pressure  dropped  off  significantly  at  larger  angles  in  this  experiment  does 
not  prove  conclusively  that  the  acoustic  wavefront  produced  by  the  laser  pulse  was 
directional. 

B.  RANGE  DEPENDENCE 

The  next  experiment  conducted  was  designed  to  test  the  expected  1/r  dependence 
of  the  acoustic  pressure  using  the  high  frequency  hydrophone.  This  experiment  was 
performed  by  moving  the  hydrophone  away  from  the  laser  spot  along  the  30°  line  where 
the  largest  pressure  response  was  achieved.  Since  a  high  speed  camera  was  not  available 
to  detennine  the  distance  between  the  laser  spot  and  the  hydrophone,  the  distance  was 
calculated  by  dividing  the  nominal  speed  of  sound  in  fresh  water  of  1480m/s  by  the 
difference  between  the  time  the  acoustic  pulse  reaches  the  hydrophone  and  the  time  of  the 
laser  pulse.  The  time  it  takes  the  laser  pulse  to  reach  the  water  is  negligible.  Therefore, 
this  approach  underestimates  the  distance  since  the  speed  of  the  initial  shock  wave  is 
higher  than  the  speed  of  a  linear  sound  wave.  However,  since  the  shock  waves  from  laser 
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pulses  of  this  magnitude  are  expected  to  dissipate  at  fairly  short  ranges,  this  should  be  a 
fairly  good  estimate  several  centimeters  from  the  laser  spot.  For  example,  if  the  shock 
wave  travels  at  4000m/s  and  dissipates  to  a  linear  wave  within  the  first  centimeter  of 
travel,  the  error  in  the  distance  estimate  would  be  about  13%  at  5cm.  The  error  decreases 
to  6%  at  10cm.  The  results  are  shown  in  Figure  15.  The  plot  shows  the  expected  roughly 
1/r  dependency. 


Pressure  vs  Range  for  High  Frequency  Hydrophone 


Range  (cm) 


Figure  15  Determination  of  the  1/r  dependence  of  the  acoustic  pulse  using  high 
frequency  hydrophone  on  a  log-log  plot. 

The  variation  of  pressure  with  range  was  also  checked  using  the  low  frequency 
hydrophone.  The  results  are  shown  in  Figure  16.  A  point  of  interest  in  this  figure  is  that 
the  pressures  are  two  orders  of  magnitude  lower  than  those  from  the  high  frequency 
hydrophone.  This  is  explained  by  the  fact  that  the  frequency  of  the  acoustic  pulse  is  well 
above  the  calibrated  range  of  the  low  frequency  hydrophone  as  seen  in  the  calibration 
curve  in  the  Appendix.  Despite  this  problem,  the  measured  pressure  still  shows  the  same 
roughly  1/r  range  dependence. 
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Pressure  vs  Range  for  Low  Frequency  Hydrophone 


Range  (cm) 


Figure  16  Determination  of  the  1/r  dependence  of  the  acoustic  pulse  using  low 
frequency  hydrophone  on  a  log-log  plot. 

Since  neither  the  spot  size  produced  from  the  laser  pulse  nor  the  size  of  the 
cavitation  bubble  could  be  determined  without  the  use  of  a  high  speed  camera,  no  direct 
comparison  can  be  made  to  the  pressures  achieved  by  other  research  groups.  Flowever,  a 
rough  comparison  can  be  made  to  Vogel’s  results  simply  based  on  the  pulse  length  and 
energy.  For  instance,  for  his  experiment  using  a  30ps  pulse  length  with  an  energy  of  lmJ 
he  got  a  pressure  of  about  lxl05Pa  at  a  distance  of  10cm.  In  the  experiments  conducted 
for  this  thesis,  using  a  pulse  length  of  about  lOOfs  and  energy  of  2mJ,  a  pressure  of  about 
2xl04Pa  was  obtained.  This  is  at  least  consistent  with  Vogel’s  work  since  he  found  the 
conversion  of  laser  energy  into  mechanical  energy  to  be  substantially  lower  for 
picosecond  pulses  as  opposed  to  nanosecond  pulses. 
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V.  MODELS 


The  modeling  effort  in  this  thesis  seeks  to  predict  the  characteristics  of  the 
pressure  wave  produced  from  the  interaction  of  a  high  energy  laser  pulse  in  water  by 
using  an  existing  finite  element  simulation  program.  For  this  preliminary  study,  only  the 
simplest  of  the  laser  parameters  are  used,  i.e.,  total  pulse  energy,  pulse  length,  and  spot 
size.  No  attempt  is  made  in  this  thesis  to  take  into  account  the  details  of  the  pulse  shape 
either  temporally  or  spatially.  It  is  simply  assumed  that  the  laser  energy  is  deposited  in  a 
spherical  volume  uniformly  over  the  time  period  of  the  pulse  length.  Likewise, 
interactions  of  the  pressure  wave  with  an  air-water  interface  are  also  ignored  at  this  point. 
The  goal  of  this  study  is  to  examine  the  output  of  the  finite  element  program  in  terms  of 
peak  pressure,  the  decrease  of  pressure  with  range,  and  the  width  of  the  pressure  pulse. 
One  of  the  biggest  problems  with  this  approach  is  the  fact  that  the  efficiency  of  the 
electro-optic  to  mechanical  energy  conversion  is  not  understood  well  enough  to  be 
modeled.  Therefore,  the  only  alternative  is  to  adjust  the  total  pulse  energy  until  the 
pressure  matches  the  experimental  value.  The  output  then  gives  a  prediction  of  the 
conversion  efficiency,  which  can  be  compared  with  experimental  values.  This  relatively 
modest  goal  shows  the  basic  ability  of  the  simulation  program  to  handle  this  extremely 
non-linear  process. 

The  first  computer  program  used  was  a  finite  element  program  called  COMSOL 
Multiphysics.  Along  with  the  basic  multiphysics  program,  COMSOL  also  has  add-on 
modules,  which  are  tailored  for  specific  applications.  For  example,  the  Acoustics 
Module  contains  built-in  application  modes  and  boundary  settings  for  the  modeling  of 
acoustic  wave  propagation  in  solids  and  fluids.  COMSOL  also  provides  the  ability  to 
modify  the  equations  of  state  of  the  materials  of  interest.  The  problem  is  that  the 
deposition  of  high  energy  laser  energy  is  highly  nonlinear.  COMSOL  could  conceivably 
handle  the  highly  nonlinear  physics,  but  it  would  have  required  heavy  modification  of  the 
physics  equations. 
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To  avoid  the  difficulty  of  having  to  input  the  equations  needed  to  handle  this  non¬ 
linear  process,  another  finite  element  program,  AUTODYN,  was  chosen  for  the  modeling 
effort.  AUTODYN  software  is  designed  for  modeling  the  nonlinear  dynamics  of  solids, 
fluids  and  gas  as  well  as  their  interactions.  Its  strong  point  is  the  modeling  of  explosions. 
The  efficient  generation  of  sound  by  laser  energy  involves  the  fonnation  of  a  plasma  with 
an  extremely  high  energy  density.  It  is  essentially  a  small  explosion.  Therefore,  it  was 
decided  that  AUTODYN  would  be  a  better  choice  for  this  application.  The  following 
sections  will  discuss  three  different  approaches,  which  were  attempted  before  finally 
settling  on  the  final  approach.  Comments  are  provided  as  to  the  difficulty  with  these 
alternative  approaches.  The  final  approach  involved  simulating  the  laser  pulse  as  a 
detonated  high  energy  explosive.  Various  results  of  this  method  are  discussed. 

A.  RADIAL  VELOCITY  APPROACH 

The  first  approach  was  an  attempt  to  model  the  results  of  Vogel  et  al.  [1996]  by 
using  the  radial  velocity  of  the  wall  of  the  cavitation  bubble.  The  6ns,  lOmJ  laser  spot 
had  a  diameter  of  7.6  pm  and  a  shape  that  was  approximately  spherical.  Vogel  predicted 
the  velocity  of  this  wall  as  a  function  of  time  by  using  the  Gilmore  model  of  cavitation 
bubble  dynamics  and  by  requiring  the  solution  to  yield  a  maximum  bubble  radius 
conforming  to  the  experimental  value.  He  also  measured  the  bubble  wall  velocity 
experimentally.  The  calculated  maximum  bubble  wall  velocity  was  1106m/s,  and  the 
experimental  value  was  approximately  twice  as  high  at  2450m/s.  This  discrepancy  was 
attributed  to  the  failure  of  the  Gilmore  model  to  model  the  shock  front  accurately. 
Interestingly,  the  calculated  shock  wave  pressure  at  1mm  was  higher  than  the 
experimental  value.  The  numerical  calculation  yielded  a  pressure  at  the  shock  front  of 
about  70MPa  while  the  measured  shock  wave  pressure  at  1mm  was  20MPa.  These 
values  agree  within  an  order  of  magnitude. 

To  begin  the  modeling  process,  a  rigid  sphere  of  the  above  mentioned  diameter 
was  created  in  AUTODYN  and  surrounded  by  water.  While  creating  this  sphere, 
AUTODYN  allows  the  user  to  select  the  equation  of  state  that  will  be  applied  to  the 
materials  in  the  problem.  In  this  case,  the  shock  equation  of  state,  as  defined  by  the 
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Rankine-Hugoniot  equations,  was  chosen  to  model  the  pressure  wave  produced  in  the 
water  as  a  result  of  this  expanding  sphere.  During  this  phase,  the  user  also  decides  how 
the  problem  will  be  meshed.  Meshing  is  the  process  of  dividing  the  materials  in  the 
problem  into  many  small  pieces  (finite  elements)  to  which  the  program  will  perform 
calculations  to  determine  the  properties  of  the  overall  problem.  For  this  problem,  a 
variable  mesh  was  used  which  divided  the  material  of  the  rigid  inner  sphere 
corresponding  to  the  expanding  bubble  and  the  water  close  to  the  sphere  into  much 
smaller  pieces  than  the  material  further  away.  This  was  done  to  minimize  the  run  time  of 
each  trial  while  maintaining  the  accuracy  of  the  calculation  close  to  the  sphere. 

AUTODYN  also  allows  the  user  to  choose  from  predefined  boundary  conditions. 
The  “Flow  Out”  boundary  condition  was  chosen  for  the  outer  surface  of  the  water  since 
this  would  most  closely  mimic  an  infinite  boundary.  AUTODYN  also  allows  the  user  to 
input  an  initial  condition  to  any  material  in  the  problem.  In  this  case,  the  material  in  the 
sphere  was  given  an  initial  condition  of  a  constant,  expanding  radial  velocity.  The  next 
step  of  the  problem  setup  is  to  put  gauges  into  the  problem  that  will  measure  the  desired 
parameters.  The  figure  produced  by  AUTODYN  for  an  initial  radial  velocity  of  1 107m/s 
is  shown  in  Figure  17  at  time  1.502  p  s  after  the  initiation  of  the  problem.  In  this  case,  a 
circular  outer  material  with  a  radius  of  2500  pm  was  created  surrounding  the  inner 
sphere  with  a  radius  of  3.8  pm.  In  the  figure  one  can  see  the  position  of  the  gauges 
(labeled  1  and  2),  the  pressure  contours,  and  the  units  used  in  calculations.  The  shock 
pressures  produced  in  this  case  were  approximately  40  kPa  measured  at  250  pm.  Using 
the  expected  1/r  relation  for  pressure  dissipation  below  lOOMPa  leads  to  about  lOkPa  at 
1mm.  These  pressures  are  approximately  three  orders  of  magnitude  lower  than  those 
measured  by  Vogel  and  his  colleagues. 

Trials  were  then  conducted  with  radial  velocities  up  to  5000m/s  with  no 
appreciable  increase  in  maximum  pressure.  The  fact  that  an  increase  in  radial  velocity 
does  not  result  in  a  corresponding  increase  in  maximum  pressure  with  this  method 
suggests  a  fundamental  flaw  in  this  approach.  Therefore,  the  decision  was  made  to 
abandon  this  method  and  try  something  different.  It  is  possible  that  this  approach  could 
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be  refined  by  using  different  equations  of  state  or  by  changing  the  radial  velocity  to 
match  the  time  dependence  reported  by  Vogel.  However,  since  Vogel  et  al.  report  that 
the  shock  wave  detaches  immediately  following  plasma  formation,  it  is  not  clear  that  the 
bubble  wall  velocity  alone  can  be  used  to  predict  shock  wave  pressures. 


AUTODYN-2D  vll.O  from  Century  Dynamics 
PRESSURE  (kPa) 


1  507e+01 

1.092e+01 

6.779e+00 

2.636e+00 

-1 .508e+00 

-5. 652 e +00 

-9.795e+00 

-1 ,394e+01 

-1.808e+O1 

-2.223e+01 

-2.637e+01 


triall 

Cycle  342 

Time  1.502E+000  ps 


Units  pm,  pg,  ps 
Axial  symmetry 


Figure  17  Example  output  from  AUTODYN 

B.  INTERNAL  ENERGY  APPROACH 

The  next  approach  was  to  investigate  the  possibility  of  adjusting  the  internal 
energy  of  the  modeled  spot  to  match  the  internal  energy  produced  by  laser  pulse.  Since 
the  units  for  internal  energy  in  AUTODYN  are  J/kg,  the  value  for  internal  energy  was 
calculated  by  dividing  the  laser  pulse  energy  by  the  mass  of  the  water  in  the  volume  of 
the  laser  spot.  The  next  step  is  to  construct  the  problem  as  described  in  the  previous 
section.  The  first  problem  encountered  here  was  in  the  choice  of  materials.  Since  only  a 
small  portion  of  the  energy  added  to  the  sphere  is  required  to  vaporize  it,  the  sphere  could 
reasonably  be  modeled  as  water  vapor.  However,  AUTODYN  does  not  have  a  pre- 
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defined  choice  for  water  vapor  in  its  materials  library.  Therefore,  the  volume  was 
modeled  with  air  since,  as  a  gas,  it  was  judged  to  be  the  best  match  for  water  vapor. 
Since  air  was  used  in  the  inner  volume,  the  ideal  gas  law  equation  of  state  was  used  for 
that  volume.  Again,  the  shock  equation  of  state  was  chosen  for  the  surrounding  water 
since  the  evolution  of  the  pressure  wave  from  this  energy  density  is  expected  to  yield 
shock  waves.  The  pressure  pulse  produced  from  this  approach  is  shown  below. 


AUTODYN-2D  vll.O  from  Century  Dynamics 


Gauge  History  (  tria!3  ) 


- (l)Gauge#  2 


Figure  1 8  Pressure  vs.  Time  for  a  gauge  500  p  m  from  center  of  laser  spot. 

This  pulse  shape  looks  promising  in  that  it  has  a  roughly  bipolar  shape.  Negative 
absolute  pressure  swings  are  realistic  when  quick  pulses  put  the  liquid  into  tension  before 
it  has  a  chance  to  react  mechanically.  However,  the  temperature  profile  for  this  run, 
shown  in  Figure  19,  is  unrealistic  in  that  it  predicts  temperatures  on  the  order  of  -107K.  It 
is  considered  likely  that  this  result  represents  a  failure  of  the  Rankine-Hugoniot  equation 
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of  state  to  model  the  full  range  of  conditions  inherent  in  this  approach.  The  physically 
nonsensical  result  might  follow  from  extrapolating  the  equation  of  state  into  regimes  in 
which  it  is  no  longer  valid. 


AUTODYN-2D  vll.O  from  Century  Dynamics 


Gauge  History  (  tria!7  ) 


Figure  19  Temperature  profile  for  internal  energy  approach. 

C.  OVERPRESSURE  APPROACH 

The  next  approach  at  the  modeling  effort  was  to  try  to  produce  an  overpressure  in 

the  laser  spot  volume.  Unfortunately,  there  was  no  direct  way  to  do  this  in  AUTODYN. 

Initial  pressure  is  not  an  input  option;  however,  density  is.  Thus,  an  attempt  was  made  to 

trick  AUTODYN  into  an  overpressure  situation  by  using  the  ideal  gas  equation  of  state  in 

an  air  bubble  with  a  high  density.  The  density  required  was  calculated  as: 

PM  n  « 

p  = - ,  Eqn38 

RT 

where  P  is  the  pressure  of  the  bubble,  M  is  the  molar  mass,  R  is  the  ideal  gas  constant, 
and  T  is  the  absolute  temperature.  The  temperature  of  the  bubble  was  calculated  by 
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assuming  that  all  the  energy  of  the  laser  pulse  had  been  transformed  into  heat  energy  and 
no  work  had  been  done  on  the  surrounding  fluid.  This  overestimates  the  temperature. 
Figure  20  shows  the  result  of  this  method  for  one  of  the  trials.  The  difference  between 
this  plot  and  the  one  obtained  with  the  internal  energy  approach  is  readily  apparent.  The 
pressure  wave  starts  with  a  negative  spike.  This  would  be  expected  from  an  imploding 
bubble,  but  it  is  not  expected  from  an  overpressure  situation.  Curiously,  however,  the 
pressure  produced  is  close  to  the  order  of  magnitude  of  the  pressure  produced  in  the 
Vogel  experiment.  This  prompted  some  further  investigation  into  this  approach. 

To  start,  pressures  between  105Pa  and  106Pa  were  converted  to  density.  The 
results  of  this  conversion  are  shown  in  the  table  below.  The  maximum  pressures 
produced  were  then  plotted  as  a  function  of  the  increasing  density.  The  results  are  shown 
in  Figure  21.  These  trials  did  not  produce  the  expected  results.  As  density  increases,  the 
maximum  pressure  of  the  resulting  wave  decreases.  These  results  imply  that  trying  to 
create  an  overpressure  in  the  laser  spot  by  adjusting  the  initial  density  actually  produces 
an  imploding  bubble.  Therefore,  this  approach  of  trying  to  trick  AUTODYN  into  an 
overpressure  situation  was  abandoned. 


Table  4  Pressure  /  Density  used  for  overpressure  trials 


Pressure  (Pa) 

Density  (g/cm3) 

1.00E+05 

7.95E-04 

2.00E+05 

1.59E-03 

3.00E+05 

2.39E-03 

4.00E+05 

3.18E-03 

5.00E+05 

3.98E-03 

6.00E+05 

4.77E-03 

7.00E+05 

5.57E-03 

8.00E+05 

6.36E-03 

9.00E+05 

7.16E-03 

1.00E+06 

7.95E-03 
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PRESSURE  (kPa ) 


AUTODYN-2D  vl  1 .0  from  Century  Dynamics 


Gauge  History  (  trial2  ) 
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-1,000 
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TIME  (|JS) 

Figure  20  Pressure  vs.  time  500  p  m  from  center  of  laser  spot  obtained  with  the 
overpressure  approach 


54 


Max  Pressure  500pm  from  Bubble  Center 


Density  of  Air  (kg/m3) 


Figure  2 1  Maximum  pressure  vs.  density  of  air  500  p  m  from  bubble  center 

D.  HIGH  ENERGY  EXPLOSIVE  APPROACH 

The  final  approach  was  an  attempt  to  use  a  high  energy  explosive  to  create  the 
overpressure  condition  desired.  This  method  consisted  of  filling  the  laser  spot  volume 
with  nitroglycerine  and  igniting  the  explosive.  Again,  a  sphere  with  a  radius  of  3.8pm 
was  used  in  an  attempt  to  model  Vogel’s  lOmJ  6ns  results.  The  resultant  pressure  pulse 
is  shown  in  Figure  22,  which  was  measured  at  37 pm,  the  distance  detennined  by  Vogel 
to  be  the  maximum  extent  of  the  plasma  radius.  Of  note  here  is  that  the  pressure  pulse  is 
no  longer  bipolar.  Rather,  the  pressure  pulse  has  the  same  shape  that  the  higher 
amplitude  pulses  showed  in  Maccabee’s  work  [1987].  At  these  amplitudes,  a  negative 
pressure  swing  cannot  occur  near  the  surface,  because  the  water  cannot  withstand  the 
tension.  Instead,  the  pressure  decreases  after  the  initial  shock  wave  passes  with  an 
exponential  tail. 
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The  first  step  of  this  approach  was  to  try  to  get  the  AUTODYN  generated 
pressure  pulses  to  be  approximately  the  same  as  those  produced  by  Vogel.  To  do  this,  the 
initial  internal  energy  of  the  explosive  was  adjusted  until  it  produced  results  that  matched 
the  measure  pressure  at  the  plasma  rim  from  Vogel’s  experiment.  As  done  before,  the 
energy  in  the  laser  spot  was  calculated  by  dividing  the  1  OmJ  laser  energy  by  the  mass  in 
the  laser  spot  volume  with  a  radius  of  3.8 pm  to  get  lxlOuJ/kg.  Then  the  initial  internal 
energy  of  the  nitroglycerine  was  adjusted  until  the  magnitude  of  the  pressure  pulse 
produced  was  approximately  the  same  as  produced  by  Vogel.  This  led  to  an  internal 
energy  of  8.7xl09  J/kg. 

The  temporal  pulse  length  predicted  by  the  model  was  also  examined.  To 
determine  the  pulse  length  from  Figure  22  the  point  where  the  pulse  decayed  to  the  1/e 
point  was  used.  From  this,  it  was  determined  that  the  pulse  length  was  about  0.5  ns 
measured  at  37  pm  from  the  center.  The  same  method  was  used  to  determine  the  pulse 
length  for  both  the  numerical  calculation  and  experimental  data  from  Vogel,  which  were 
about  15ns  and  200ns  respectively.  In  Vogel’s  work,  the  numerical  calculation  was  for 
the  pressure  inside  the  bubble  during  the  initial  phase  of  bubble  expansion,  and  the 
experimental  data  was  measured  10mm  from  the  emission  center.  It  is  expected  that  a 
pressure  pulse  will  spread  out  as  it  travels,  but  measurements  were  not  taken  that  far  from 
the  emission  center  in  this  model.  Further  work  is  required  to  examine  the  model 
parameters  which  contribute  to  the  pulse  length  and  to  determine  the  expected  pulse 
length  at  greater  distances  from  the  center. 
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AUT0DYN-2D  vll.O  from  Century  Dynamics 


Gauge  History  ( tria!7  ) 


-  (l)Gauge#  2 


Figure  22  Pressure  pulse  from  high  energy  explosive  trial  for  the  lOmJ  6ns  case  for  a 
sphere  with  a  radius  of  3.8  p m  measured  at  37  p m  from  the  center  (edge  of 
the  plasma  rim  determined  by  Vogel). 

To  get  the  efficiency  of  this  model,  the  internal  energy  of  the  nitroglycerine 
required  to  achieve  the  desired  pressure  magnitude  was  divided  by  the  energy  of  the  laser 
spot,  which  resulted  in  the  same  pressure.  The  result  was  an  efficiency  of  9%.  This 
efficiency  is  comparable  to  the  6-13%  efficiency  from  Vogel  and  the  10%  efficiency 
from  Jones.  As  one  can  see  by  comparing  Figure  23  and  Figure  24  the  pressures  from  the 
model  and  from  Vogel  are  both  on  the  same  order.  Figure  23  also  shows  the  roughly  1/r" 
dependence  that  is  expected  of  shock  wave  dissipation  at  pressures  greater  than  lOOMPa. 
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Pressure  vs  Range  for  Internal  Energy  of  SJ^IO^J/kg 


2 

Figure  23  Pressure  vs.  range  from  AUTODYN  model  showing  1/r  behavior  from 
shock  wave  dissipation. 


Figure  24  Pressure  vs.  Range  from  Vogel  for  6ns  pulse  of  lmJ  and  lOmJ. 
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VI.  CONCLUSIONS 


In  the  experimental  portion  of  this  thesis,  a  couple  of  fundamental  characteristics 
of  the  acoustic  source  were  investigated.  First,  the  directionality  of  the  source  was 
considered.  At  first  glance,  it  seemed  that  the  source  was  quite  directional.  Further 
investigation  led  to  the  conclusion  that  the  reason  for  the  directionality  of  the  source  was 
most  likely  due  to  interference  between  the  direct  acoustic  pulse  and  its  reflection  off  the 
surface. 

The  next  result  of  the  experimental  section  was  the  determination  of  the  range 
dependence  of  the  acoustic  source.  Five  measurements  of  the  peak  pressure  amplitude 
were  taken  at  ranges  varying  from  about  3  cm  to  24  cm  using  both  high  and  low 
frequency  hydrophones.  The  dependence  of  the  peak  pressure  amplitude  with  range  was 
measured  as  r  12  with  the  high  frequency  hydrophone  and  as  r-0'94  with  the  low 
frequency  hydrophone.  These  measurements  showed  that  the  laser  generated  acoustic 
source  produced  the  approximately  1/r  range  dependence  expected  for  pressures  under 
lOOMPa. 

The  modeling  aspect  of  this  thesis  also  produced  some  interesting  conclusions. 
To  start,  it  was  determined  that  AUTODYN  handled  the  nonlinear  aspects  of  this  thesis 
better  and  was  easier  to  use  than  COMSOL.  It  also  showed  that,  although  there  were 
many  possible  approaches  within  AUTODYN  to  model  the  acoustic  source,  it  was 
determined  that  the  best  results  came  from  modeling  the  laser  generated  acoustic  source 
as  a  sphere  of  nitroglycerine  with  adjusted  internal  energy.  This  method  led  to  a  couple 
of  positive  results.  The  first  was  the  calculated  efficiency  was  comparable  to  that  of  the 
main  references  used  in  this  thesis.  Comparing  the  adjusted  internal  energy  to  the 
calculated  internal  energy  of  the  laser  spot  led  to  an  efficiency  of  8.7%.  This  efficiency  is 
comparable  to  the  6-13%  efficiency  from  Vogel  and  the  10%  efficiency  from  Jones.  The 
model  also  showed  the  range  dependence  of  the  peak  pressure  amplitude  of  the  pulse  to 
be  r  2  2  from  10  to  100  pm.  This  is  in  agreement  with  the  roughly  1/r2  relationship  of 
pressure  vs.  range  that  is  expected  in  the  strong  shock  wave  regime.  This  method  also 
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produced  a  pulse  length  of  about  0.5ns  measured  at  37  pm,  which  is  the  maximum  extent 
of  the  plasma  as  detennined  by  Vogel.  This  value  for  pulse  length  is  within  about  an 
order  of  magnitude  of  the  15ns  pulse  length  calculated  inside  the  plasma  produced  by  the 
numerical  analysis  done  by  Vogel. 
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VII.  RECOMMENDATIONS  FOR  FOLLOW-ON  WORK 


Although  this  thesis  shows  the  promise  of  a  laser  generated  acoustic  source  there 
is  still  a  lot  of  work  that  needs  to  be  done.  One  of  the  questions  that  arose  from  the 
theory  section  had  to  do  with  the  plasma  frequency.  If  the  energy  density  is  too  high, 
electron  densities  may  be  high  enough  to  exceed  the  plasma  frequency.  In  this  case,  all 
remaining  electro-optic  energy  would  be  reflected.  Thus,  the  upper  limit  to  the  energy 
density  is  determined  by  the  point  where  the  plasma  frequency  is  reached.  Further 
investigation  needs  to  be  done  to  determine  the  highest  energy  density  before  plasma 
frequency  is  reached,  and  how  does  the  pulse  shape  affect  the  point  where  the  plasma 
frequency  is  reached? 

Another  point  for  further  investigation  will  be  the  interaction  between  multiple 
pulses.  Due  to  the  need  for  an  acoustic  signal  in  the  hundreds  of  kilohertz,  work  will 
need  to  be  done  to  detennine  how  each  successive  laser  pulse  will  interact  with  the 
bubbles  and  heating  produced  by  previous  laser  pulses.  This  effort  will  need  to  include  a 
study  of  the  refraction  of  sound  in  the  water  heated  from  successive  laser  pulses  and  also 
the  potential  scattering  of  acoustic  waves  off  of  the  bubbles  formed. 

The  next  area  for  follow-on  work  comes  from  the  modeling  aspect.  Work  needs 
to  be  done  to  detennine  if  modeling  the  laser  generated  acoustic  source  as  a  sphere  fdled 
with  an  explosive  is  the  most  realistic  means  of  modeling  this  problem.  The  model  also 
needs  to  be  expanded  so  that  the  pressure  pulse  can  be  measured  much  further  from  the 
laser  spot  to  determine  if  the  pulse  length  result  produced  by  Vogel  at  10mm  can  be 
reproduced.  Once  that  is  detennined,  the  next  step  will  be  to  determine  a  way  to  model 
multiple  pulses  and  to  see  if  AUTODYN  can  predict  the  way  multiple  pulses  will  interact 
with  each  other. 
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Finally,  the  question  of  efficiency  must  be  addressed.  For  example,  the  efficiency 
of  the  free  electron  laser  is  expected  to  reach  about  10%.  Combined  with  a  conversion 
efficiency  for  the  electro-optic  to  acoustic  process  of  approximately  10%,  this  leads  to  an 
overall  efficiency  of  1%.  Further  studies  of  the  potential  advantages  of  a  laser  source  are 
required  to  determine  whether  they  outweigh  the  disadvantage  of  this  low  efficiency  and 
to  compare  the  cost  with  more  traditional  acoustic  sources. 
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APPENDIX 


HYDROPHONE  SENSITIVITY 


TC4034-1 
3006287 
TC4034 
2006-10-28 
4492,  49 

PHO  0  250Hz:  -217.0dB. 


AmpStude:  10.0  Vrms 

Pulse  Width:  428.6  ps 

Rep  Rate:  33.3  ms 

Averages:  8 


Temperature: 
Depth: 
Distance: 
Tested  by: 


21.23*0 
15  m 
0.50  m 

HHP*a V*toCf 


g  7  g  «  io 


■M  M  HO  "  Stt 


Figure  25  Low  Frequency  Hydrophone  Sensitivity.  Sensitivity  in  dB  re  lV/pPa  vs. 
frequency  on  a  linear-log  plot. 


63 


Frequency  response  of  hydrophone  system  comprising 
Prea  mpHfier  SN:  W27B003,  DC  Coupler  SN:  272  and  1.0  mm 
Needle  Hydrophone  SN:  813 


Time  domain  waveform  in  non-linear  field  of  hydrophone  system 
comprising  Preamplifier  SN:  W 279003,  DC  Coupler  SN:  272  and  1.0 
mm  Needle  Hydrophone  SN:  813 

800  - 


-400  4 - T - . - . - . - » - - - \ 

288  290  292  294  296  298  300  302 


Time  (us) 


Figure  26  Sensitivity  of  Needle  Hydrophone  with  serial  number  813. 
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Table  5  Calibration  data  showing  sensitivity  and  uncertainty  at  each  frequency  for 
hydrophone  with  serial  number  813. 


FREQUENCY  RESPONSE  CALIBRATION  DATA 


Device  Identification 


In-House  Preamplifier  and  DC  Coupler 


Preamplifier  SN 


DC  Coupler  SN 


W279003 

272 


Hydrophone  Diameter 
Hydrophone  Type 
Hydrophone  SN 


1.0  mm 
Needle 
813 


Frequency  (MHz) 

Sensitivity  (mV/MPa) 

Uncertainty(%) 

1 

1225.6 

14 

2 

819.1 

14 

3 

878.8 

14 

4 

916.2 

14 

5 

908.6 

14 

6 

903 

14 

7 

687 

14 

8 

892.3 

14 

9 

686 

15 

10 

887.2 

15 

11 

865.2 

16 

12 

865.6 

17 

13 

817.3 

18 

14 

739,2 

19 

15 

672.2 

20 

16 

582.3 

22 

17 

498.3 

25 

18 

428.9 

32 

19 

335.3 

34 

20 

203 

38 
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Table  6  Calibration  data  showing  sensitivity  and  uncertainty  at  each  frequency  for 
hydrophone  with  serial  number  825. 


FREQUENCY  RESPONSE  CALIBRATION  DATA 


Device  Identification 


In-House  Preamplifier  and  DC  Coupler 


Preamplifier  SN  W279003  Hydrophone  Diameter  1.0  mm 

DC  Coupler  SN  272  Hydrophone  Type  Needle 

Hydrophone  SN  825 


Frequency  (MHz) 

Sensitivity  (mV/MPa) 

Uncertainty(%) 

1 

1441.6 

14 

2 

1081 

14 

3 

1094.1 

14 

4 

1006.8 

14 

5 

1050.2 

14 

6 

1011.7 

14 

7 

1027.6 

14 

8 

1042.8 

14 

9 

998.5 

15 

10 

1088.1 

14 

11 

1092.7 

15 

12 

1083.5 

15 

13 

1057.9 

15 

14 

1005.9 

15 

15 

954.6 

15 

16 

887 

16 

17 

827.2 

18 

18 

747.7 

18 

19 

654.2 

19 

20 

592 

19 
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Sample  Directivity  for  1.0mm  Needle  Hydrophone 


Angle  (degrees) 


Figure  28  Directivity  for  needle  hydrophones 
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Gain  in  dB 
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Gain  lit  dB 
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